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COMPUTER PROGRAM GRADE II 
FOR THE DESIGN AND ANALYSIS OF HEAT-PIPE WICKS 

1.0 INTRODUCTION 

This user's manual describes the revised version of the computer 

( 1 ) 

program GRADE , which designs and analyzes heat pipes with graded- 
porosity fibrous slab wicks. The revisions, which are based on work done 
under contract NAS 2-8310 with NASA Ames Research Center and reported in 
Reference (2), were incorporated so that the mathematical model more com- 
pletely describes an actual graded-porosi ty-wi ck heat pipe. In particular, 
GRADE II now Includes: • - 

• Automatic calculation of the minimum condenser-end stress that 
will not result in an excess-liquid puddle or a liquid slug in 
the vapor space, 

• Numerical solution of the equations describing flow in the cir- 
cumferential grooves to assess the burnout criterion, 

• Calculation of the contribution of excess liquid in fillets and 
puddles to the heat- transport, 

• Calculation of the effect of partial saturation on the wick per- 
formance, 

• Calculation of the effect of vapor flow, which includes viscous- 
inertial interactions. 

In addition to these extended capabilities, the new version retains the 
capabilities of the original program: 

• Calculation of the optimum porosity variation and the correspond- 
ing maximum heat- transport raite, 

• Calculation of the maximum heat- transport rate at other than the 
wick's design condition (different temperature, elevations, gravi- 
tational field, etc.), 

• Calculation of the maximum heat-transport rate for a specified 
porosity distribution, which includes a uniform-porosity wick. 
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• The heat pipe can have multiple sections each having different 
tilts, 

• Multiple heat input and output zones, 

• Calculation of the total fluid charge. 

The theoretical basis for GRADE II is described in Section 2.0 and 
the instructions for preparing the input are given in Section 3.0. If 
excess-liquid effects are to be included in the calculations, a separate 
program FILLET must be run, whose output is a binary file that becomes 
additional input for GRADE II. This program is described in Section 4.0. 
Two sample programs are described in Section 5.0. The Appendix contains 
descriptions and listings of GRADE II and FILLET. 

2.0 THEORETICAL BASIS FOR GRADE IT 

This section describes the application of fundamental theoretical 
and experimental results for caoillary flow through fibrous media previ- 
ously reported in Reference (3) and (4) to the optimum design of heat- 
pipe wicks. To be specific, we consider a heat pipe, as depicted in 
Figure 1, with a fibrous slab wick for the axial transport of liquid and 
circumferential grooves for the transport across the evaporation and con- 
densation surfaces. The slab wick is ultimately limited in he at- transport 
capacity because the factors that affect its performance, the capillary- 
pressure limit and the permeability, are related inversely. Any change 
in the wick structure that increases Its capillary-pressure limit de- 
creases its permeability and vice versa. The simple uniform-porosity 
wick is optimized by selecting the fiber diameter and porosity that 
maximizes its heat transport. Such a wick, however, has an unnecessarily 
low permeability, everywhere along its length except where it begins to 
dry out under maximum load. A further capacity increase is possible if 
one considers a wick whose porosity varies along its length. With a 
graded-porosi ty wick, the porosity *s optimally varied such that at 
every axial location *t is only as low as required to ensure the wick 
remains nearly saturated. Thus, the permeability is everywhere as high 
as possible. The potential increase in capacity over a uniform porosity 
wick depends on the particular applications, but it Is often greater than 
a factor of two. 


2 
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Figure 1. Cross-Section of a Fibrous-Slab-Wick Heat Pipe 5 

-3- ; 



I 


i 


26263-6026-HU-00 


V4. 

V 

\ 


t 





/ 

9 

/ 


1 


4 


f 

i 


2.1 PROPERTIES OF A FIBROUS MICK 


Me first summarize the results of Reference (3) and (4) for capillary 
flow through porous media and then derive the relationships that are re- 
quired for the design of a wick. Expressions for the capillary-pressure 
limit P c and the permeability K Q for a wick of uniform porosity < consist- 
ing of a three-dimensional random distribution of fibers of diameter 5 
were shown to be 

P c = 3.2465 H(o/6)(l - c)/t (1) 


and 

Kq * (3/8)6 2 [e/(! - e)]/ 

< - h 

X 4(1 - e ) - (1 - tr - 2 1 n( 1 - c) - 3 

' 10(7 - 0 + [1 - 0 - e) 2 ]/[l + (1 - e) 2 ]J 


( 2 ) 


where a is the surface tension and H is the hysteresis constant that is 
unity if the liquid front is advancing in the wick and an empirically 
found value of 1.955 if it is receding. Actually, the wick does not 
empty abruptly when the capillary-pressure limit is exceeded, but rather 
it progressively desaturates. The wick is envisioned as consisting of 
local regions having porosities that are normally distributed with a mean 
value c Q and a standard deviation o^. The fraction of the wick with a 
porosity that lies between c and e ♦ de is given by 


f( 




7 d) 


- (< 


to) 2 / 2 * 7 / 


.Tr 


(3) 


°d 


The standard deviation was found experimentally to correlate with the 
mean porosity by the expression 

o d * 0.22(1 - e 0 ) (41 

When the wick is subject to a vapor-liquid pressure difference P» 
which we call the capillary stress, a local region Is filled with liquid 
if its porosity is sufficiently low that the capillary-pressure limit 
given by £q. (1) exceeds the capillary stress. The saturation fraction. 
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which is the ratio of the liquid content of the wick to the content when 
it is completely saturated, was shown to be 

S - F[e* - e 0 )/c d ] - (o d /e 0 )f[(e* - c Q )/a d ] (5) 

where f(z) is the standardized normal distribution and f(z) is the stan- 
dardized cumulative distribution, and e* is the critical value of the 
local porosity for which the capillary pressure limit equals the capillary 
stress. Its value, obtained from Eq. (1), is 

P* /r* “1 r«i 


«* - « ♦ (6) 

To obtain an expression for the permeability of the partially saturated 
wick, Eq. (2) is applied to those regions with a porosity below the 
critical value. The resulting expression is 


K(c 0 , e*» 6, o d ) * £ K 0 ( 6 » e)f(e, eo» c d )de (7) 

To this point, we have sunmari zed the results of Reference (3) and 
(4). We now calculate the mean porosity that maximizes the permeability 
for a prescribed capillary stress. If the porosity is too high, the 
wick is unable to hold liquid at the prescribed stress which results in 
a low permeability. If, on the other hand, the porosity is too low, the 
wick will remain nearly saturated, bit the fibers are unnecessarily close 
together which also results in a low permeability. Equation (7) is the 
basis for the optimization. A specified value of the capillary stress 
fixes e* by way of Eq. (6). The dependence of the permeability on the 
fiber diameter & is eliminated from Eq. (7) by using 5^ to nondimension- 
a'ize K, and the dependence on S d is eliminated with Eq. (4). The result- 
ing expression for the dimensionless permeability K = K/6 2 is 


K(e„, 


4{1 - c)-(l - e) 2 - 2 ln(l - c) 
10.3(e - cj 2 


ln(l - c) + 


[1 - e 


1 +0 
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Figure 3. Key Relationships for Design of a 
Graded- Porosity Hick 
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This Integral was evaluated numerically with Simpson's rule on a computer. 
The results are shown in Figure 2. Also shown are lines of constant 
saturation fraction, which were calculated from £q. (5) witf given by 
Eq. (4), and the locus of points for which the permeability is maximum. 

Figure 3 relates the critical porosity .* to the dimensionless cao- 
illary stress P :/; and displays the optimum mean porosity ind the cor- 
responding saturation fraction as a function of the stress. A graded- 
porosity wick ideally would have a porosity variation that follows the 
optimum-porosity curve as the capillary stress builds fro>« a low value 
at the end of the condenser to a high value at the end of the evaporator. 
The optimum saturation fraction would range from a value above 0.8 at a 
low stress typical of the condenser region to a value below 0.7 at a high 
stress typical of the evaporator region. The fact that such a wick oper- 
ates with a liquid fill well below that required to saturate it presents 
practical problems. If, for example, a fluid charge is used that is 
sufficient to completely saturate the wick, then at the maximum heat- 
transport rate liquid will be given up that could result in flooding of 
the condenser. If, on the other hand, a fluid charge is used that is 
just, sufficient to provide the optimum saturation fraction at the maximum 
heat- transport rate, then there is no guarantee that the liquid will be 

properly distributed along the wick. To avoid these problems, the wicks 

are designed with a porosity variation that provides a uniform high-level 
of saturation. Thus, instead of operating along the peaks of the parti- 
ally saturated permeability curves of Figure 2, the wick is designed to 
operate to the left of the peaks along a line of constant saturation 
fraction. Equation (5) is usbd to obtain the expression for the porosity 
that provides the desired saturation fraction. The equation is trans- 
cendental in £ 0 , and it must be solved iteratively. For hi ah levels «f 

saturation, however, the second term of Eq. (5) is small compared to the 

first, and an accurate approximation for *o can be obtained by neglecting 
It, which results in 

c 0 . - e:'<s?/4.s . ( 9) 

1 r F \$}/4.5 

where we have used Eq. (4) to eliminate o^. Equation (9) was used to 
calculate the curve of Figure 3 which gives the porosity as a function of 
capillary stress that provides a saturation fraction of 0.9. [The criti- 
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cal porosity is first calculated as a function of stress fro® £q. (6)3. 
2.2 DESIGH OF A WICK 


We nov focus attention on the hydrodynamic optimization of a heat 
pipe for maximum heat transport in a given duplication. One must be 
alert, however, to tie possibility that other limiting factors may come 
Into play before the hydrodynamic wicking limit Is reached. Such factors 
are, for example, the heat-flux limit due to boiling in the wick, and 
the sonic vapor-flow limit. The procedure described herein is used to 
calculate the optimum wick porosity variation for a fixed beat-pioe geo- 
metry. For Vie heat-pipe diameter considered, a change in wick area may 
further Increase the capacity. If the capillary stress is due primarily 
to liquid flew through the wick, an increase in wick area will increase 
the capacity. If, on the other hand, the stress Is due primarily to 
vapor-flow pressure drop, or if there is a relatively large vapor-space 
capillary back pressure, which we will see presently can adversely affect 
the porosity variation, a reduction in wick area will increase the capa- 
city. In fact, for a given heat-pipe diameter, there is always an opti- 
mum wick area. 


The key equation describing the heat-pipe hydrodynamics governs the 
axial variation of the capillary stress. In a gravitational field, how- 
ever, it varies hydrostatically across the heat pipe as well; so to have 
a unique value at every axial location, we take its value at the top of 
the wick. The equation governing the rate of Increase of stress P with 
axial distance x from the condenser end Is 
jn v, m{x) .. 

& ’ ” Cv) 9 <& * 

~2 

7 PyU J 

~2d" ' F * d? (10) 

The first term on the right of the equal sign gives thre stress increase 
due to liquid of kinematic viscosity vj flowinq through a wick of cross- 
sectional area A* at a mass rate m{x). The permeability K depends on 
the local porosity c 0 ««** on the capillary stress through the critical 
porosity c*. 


• 10 - 
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The second term gives the change in hydrostatic pressure in the 
liquid of density due to changes of heat-pipe elevation h(x) measured 
from a horizontal reference plane to the top of the wick. 

The third and fourth terms are due to vapor flow, as discussed in 
Reference (2). Here IT is the average velocity in the vapor space. Re is 
the Reynolds number, and D is the hydraulic diameter of the vapor space. 
The third term is due to viscous shear on the walls, and the fourth is 
due to inertial effects. K and F s are, respectively, an average friction- 
factor coefficient and shape factor, which are calculated according to 
Reference (2) to give the proper balance between inertial and viscous 
effects. The mass flow rate is related to the latent heat of vaporiza- 
tion hfg and the heat input per unit length Q(x) (assumed negative in 
regions of condensation) by 

m(x) = - (1/hfg) ' Q(x) dx (11) 

The ‘optimum porosity distribution is calcualted by numerically inte- 
grating Eq. (1) with an assumed value for the heat load. At each step of 
the integration, the critical porosity e* and the wick porosity e 0 are 
calculated from Eqs. (6) and (9) with a specified high level of satura- 
tion fraction S and with the hysteresis constant H = 1.955 for liquid on 
the verge of emptying. Because of hysteresis, however, the calculated 
porosity s 0 may be too high for the wick to self-fill to the specified 
level of saturation under a zero heat load. Therefore, Eqs. (6) and (9) 
are used again with H = 1 for liquid filling the wick and the stress 
given by integration of Eq. (10) with m(x) = 0 to calculate c* and e 0 
for the wick to fill. These latter values are used if the porosity 
required tc fill the wick under zero load is lower than the porosity 
required to sustain the stress under the assumed load. 

Once e 0 and c* have been determined, the value of the permeability 
at the particular integration step is calculated from Eq. (8). In 
regions of evaporation, the circumferential grooves are checked at each 
step to see whether or not they dry up. The subroutine DRY, which is 
based on the mathematical model of Reference (2), is called to make this 
check. It takes into account viscous flow in the groove under the action 
of surface tension and gravity. If the grooves are found to dry up. 


- 11 - 
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the integration is stopped, the assumed heat load is reduced, and the 
integration is repeated. If the integration continues to the evaporator 
end of the heat pipe without groove dry-up occurring, the assumed heat 
load is increased, and the integration is repeated. A binary search is 
used to find the maximum heat load that does not result in groove dry-up. 

The calculation of the condenser-end stress used to begin the inte- 
gration is crucial. The stress must be high enough to prevent a liquid 
puddle or slug from forming in the lowest vapor space. If, however, the 
stress is set too high, the wick must begin with an unnecessarily low 
porosity to enable the wick to fill. In the condenser region, where the 
wick porosity is relat’vely high, a small reduction in porosity jean 
result in a large reduction in permeability. For example, if a satura- 
tion fraction of 0.9 is specified, then we see from Figure 2 that a 1 
percent reduction in a typical condenser-end porosity of 0.88 leads to 
a 15 percent reduction in the permeability. Therefore, for a high heat- 
transport capacity, the condenser-end stress should be kept as low as 
possible. This is one reason why, as discussed previously, a reduction 
in wick area can result in an increase in capacity. The increased size 
of the vapor space reduces its capillary back pressure and allows a 
higher condenser-end porosity. 

The first requirement on the condenser-end stress is that it. must be 
high enough that a puddle does not form in the lower vapor space. For 
simplicity, we restrict our attention to the situation depicted in Figure 
1 where the slab wick is horizontal. If the capillary stress at the top 
of the wick has been increased to a point where a puddle is just about to 
disappear, the radius of curvature of the meniscus of the puddle is nearly 
equal to the tube radius R. The stress c/R at the puddle surface is 
related hydrostatically to the stress at the top of the wick; thus the 
stress required to prevent a puddle is 

P 0 * o/R + (p a - p v ) g h w (12) 

where h w is the distance between the top of the wick and the bottom of 
the tube. The second requirement on the condenser-end stress is that it 
should be high enough to prevent a liquid slug in the lower vapor space. 

- 12 - 
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The conditions under which a slug will form, presented in Reference (2), 
are calculated by subroutine VSBKS. The stress used to begin the integra- 
tion is the greater of the value for the formation of a puddle and a 
liquid slug. 

3.0 INPUT FOR GRADE II 

The input is in Fortran NAMELIST form. The required parameters are 
defined and discussed below. An input form is given in Table 1. 

3.1 HEADINGS 

After writing on the first card or line the NAMELIST identifier 
$GRDATA, the user then inputs two lines of descriptive information by 
writing on one card or line HD1 = 60H followed by up to 60 characters of 
title and on the next card or line HD2 ■ 60H followed by another 60 
characters. GRADE II will print these two lines at the beginning of the 
output. 

3.2 FLUID PROPERTIES 

GRADE II automatically computes the required fluid properties for 
one of several fluids, which the user specifies by selecting a value of 
LIQ from the following list: 



Fluid 

Temperature Range 

LIQ ■ 1 

Water 

(OC < T < 204C) 

LIQ = 2 

Ammonia 

(-78C < T < 88C) 

LIQ = 3 

Methyl Alcohol 

(-96C < T < 1930 

LIQ = 4 

FREON-21 

(-48C < T < 152C) 

LIQ * 5 

Ethane 

(-93C < T < 27C) 

LIQ = 6 

Methane 

(-173C < T < -84C) 

LIQ = 7 

Nitrogen 

(-207C < T < -157C 


The properties are for a temperature TKELYN that the user inputs in 
degrees Kelvin. All other fluid properties are automatically computed 
for that temperature. If another fluid is used, set LIQ = 0. Then, 
y values must be specified for the following quantities: 
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Quantity 

Symbol 

Units 

Liquid density 

RHOL 

Kg/cu. m 

Vapor density 

RHOV 

Kg/cu. m 

Liquid viscosity 

VI SL 

N’S/sq. m 

Vapor viscosity 

VISV 

H*S/sq. m 

Surface tension 

ST 

N/m 

Latent heat 

HFG 

J/Kg 


3.3 GEOMETRICAL PARA.METERS OF THE HEAT-PIPf CROSS SECTION 

As shown in Figure 1, the heat pipe uses a slab wick In either a 
horizontal or vertical orientation. The input parameters that specify 
the cross-sectional geometry are: 


Cuanti ty 

Symbol 

Uni ts 

Tube inside diameter 

HPID 

cm 

Slab-wick thickness 

WKTH 

cm 

Geometrical parameter: 
Horizontal slab wick 

IGEOM » 0 


Vertical slab wick 

IGEOM » 1 



All other parameters of the heat-pipe cross section, such as the wick 
area, vapor-space hydraulic diameter, etc. are automatically calculated. 

3.4 HEAT INPUT 

The user specifies the heat-input distribution by specifying the 
fraction of the total heat-transfer rate for up to ten segments of the 
heat pipe. Heat transfer is assumed to be constant along the given seg- 
ment. Values for the following parameters are required: 

NQ - The number of segments, which must not exceed ten, into 
which the heat pipe is divided. 

XQ(I) - The length of the I th segment in cm. The segments must 

be nunbered consecutively along the heat pipe beginning at 
the condenser end. 

FQ(I) - The fraction of the total heat- transfer rate entering the 
Ith segment. If the I th segment is a condenser, FQ(J) is 
negative. If the I th segment is adiabatic, FQ(I) is zero. 


% 


- 16 - 
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QD0T - A nominal heat-transfer rate In watts* which is the user's 
best guess at the maximum. A close guess reduces the number 
of iterations to the final answer. 

3.5 ELEVATIONS 


The user specifies the heat-pipe orientation in a gravitational field 
by Inputing values for elevations of points along the heat pipe where the 
slope changes. Between points. GRADE assumes a linear variation of eleva- 
tion. Values of the following parameters are required (except for zero 
gravity): 

NEIEV - The total number of points along the heat pipe, which 
must not exceed 10. 

XELEV(l) The distance along the heat pipe In cm to the I th point. 

Both ends of the heat pipe must be input, therefore the 
first point must be at zero distance [XELEV(l) * 0.0], 
and the last at the total heat-pipe length CXELEV(NELEV) « 


LJ. 

ELEV(I) - The elevation In cm of the point relative to a hori- 
zontal reference plane. 

GEE - Gravitational acceleration in standard gravities. 

3.6 CIRCUMFERENTIAL GROOVE PARAMETERS 


The cross-section of the circumferential grooves is trapezoidal as 
shown in Figure 4. 
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The following parameters must be input: 


Quantity 

Symbol 

Units 

Number of grooves per cm 

GRVS 

cm'l 

Groove opening 

W 

cm 

Groove depth 

DEPTH 

cm 

Groove half angle 

PHI 

degrees 

Wetting angle 

AMGWET 

degrees 


3.7 WICK PARAMETERS 

The user can use the program either to design an optimum graded- 
porosity wick and compute its capacity, or he can use it to compute the 
capacity of a wick with a specified porosity distribution. i 

3.7.1 De sign of a Graded-Porosity Wick 

To design a graded-porosity wick, the user must specify: ! 

Quantl tv Symbol 

Fiber diameter DIAF 

Saturation fraction S 

Minimum allowable porosity EPSMIH 

The saturation fractions are the uniform high level of saturation 
that the wick is to maintain. We have been using S * 0.9. When the 
porosity variation Is to be designed, the parameters LASTEPS and HEPS 
must be set to zero or, equivalently, just not included in the NAMELIST. 
The mi pi mum porosity EPSMIH is set so that a wick will not be designed 
that is too dense to manufacture. 

3.7.2 Capacity at Off-Ootlmum Operation 

The program is set up to run several cases: one NAMELIST input is 
simply followed by another. When a wick is designed and the user desires 
to calculate the performance of that wick at off-optimum conditions (a 
different temperature or evaporator elevation, for example), he Sicily 
sets the parameter LASTEPS ■ 1. This causes the porosity variation to 
be that of the previous case (be sure to set HEPS » 0) . 


Units 

cm 




1 


I 




l 


% 


1 

I 
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3.7.3 Specified Porosity Distribution 

If the user desires to confute the capacity of a wick with a speci- 
fied porosity distribution, he specifies the local porosity of a number 
of points along its length. Between points, values of the porosity are 
calculated by linear interpolation. The required Input parameters are: 


Quantity 

Symbol 

Units 

Number of porosity points (up to 10) 

.HEPS 


Distance to I th point 

XEPS(I) 

cm 

Porosity of I th point 

EPSX(I) 


Hick fiber diameter 

DIAF 

cm 


The first point must be at the condenser end fXEPS(l) * C] and the last 
point must be at the evaporator end UEPS(HPHI) ■ heat-pipe length]. 

Set $ • 0 and IASTEPS » 0. 




* 



t 


s 


s 



i 



3.8 OTHER INPUT PARAICTERS 

OX - The integration step size in cm. 

IPRIMED - Equals 0 if the user requires the wick to self-prime 
under no load at the operating elevations. 

- Equals 1 if the wick is allowed to self-prime level under 
no lead before the heat pipe is raised to the operating 
elevations. 

ROUGH - The average surface roughness of the vapor spaces, which 
is used for the calculation of the turbulent friction- 
factor coefficient. 

I FITS * Equals 1 if the contribution of excess-liquid fillets 

and puddles are to be included, in which case an additional 
Input file is needed (see Section 4). 

- Equals 0 If the excess-liquid contribution is not to be 
Included. 

NCASE - Equals 1 if a HAMEL! 5 T input for another case is to follow, 
which is exactly like the first except only those parame- 
ters that are to be different in the new case are included. 

- Equals 0 if the present case is the last case. 

SEND - Ends present HAMEL 1ST input. a 

■ ; : - ' . ' ■ ■ ’■ .. ' _ •' 
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4.0 PIE PROGRAM FILLET 

If the user includes the effect of excess-liquid fillets and 
puddles, he sets 1FLTS > 1 in the HAMEL I ST input for GRADE II. Ihis 
causes GRADE II to read data from a binary file. TAPE 7, which is the 
output of the program FILLET. The theoretical basis for FILLET is de- 
scribed in Reference (2). In essence, it nuaerically integrates the 
differential equations that describe the free -surface shape of fillets 
and puddles that can exist in the heat pipe. The output fro* - . FILLET is 
a table of total cross-sectional area and hydraulic diameter of the ex- 
cess liquid as a function of stress. t 

4.1 INPUT TO FILLET 

FILLET also uses NAMELIST input. An input form is given in Tpble 
II. The NAMELIST identifier is SFILLETD. The input variable, which have 
the sane definitions as those for CRADE II are: 


LIQ 

TKELVN 

RHO* 

ST* 

IGEOM 

WKPi 


Refer to Section 3.2 


Refer to Section 3.3 


GEE Refer to Section 3.5 

•RHO, which is the difference between the liquid and 
and vapor density (RHOL-RHOV), and the surface tension 
ST are input only if LIQ * 0. 

The NAMELIST ends with the line SEND. 
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5.0 SAMPLE CALCULATIONS 

In this section we describe two sample problems. The first is 
selected to illustrate the option for including excess liquid, while 
the second is selected to illustrate the design of a graded-porosity 
wick. 

5.1 A SIMPLE METAL-FELT SLAD-WICk HEAT PIPE 

The input file for the first heat pipe we are considering is given 
in Table III. The heat pipe is 30-ca long, with condenser, adiabatic 
and evaporator lengths of 10 cm each (XQ • 3* 10.). The inside diameter 
of the tube is 1. cm (HPIO ■ 1.), the wick thickness is .4 cn (WKTH * .4) 
and it is vertical (IGEQM ■ 1). The wick is a slab of felt metal which 
has a fiber diameter of 0.002 cm (OIAF • .002). To specify a uniform 
porosity, we specify the porosity of two points (HEPS ■ 2) at the con- 
denser end and the evaporator end (XEPS * 0., 30.), At ?ach point we 
set the porosity to 0.80 (EPSX • 2* .80), and thus a linear interpolation 
for points in between results in the desired uniform porosity. The cir- 
cumferential grooves have a width at the top of 0.015 c® (H * .015), a 
depth of 0.015 cm (DEPTH » .015) and a half angle of 20 degrees (PHI • 
20.). We specified 100 grooves per cm (6RVS » 100.). Upon reflection 
one sees that it would be impossible to cut such grooves so close to- 
gether; however , they will serve for purposes of this illustration. 
Several cases are run with the evaporator elevated 2., 0., 4., and o. cm 
higher than the condenser end. The fluid is ammonia (UQ * 2) at 30QX 
(TXELVH * 300.). 

Since the effect of liquid fillets are to be included. 1FLTS is set 
equal to unity, and the program FILLET must be run first. The NAMELIST 
input to FILLET is shown in Table IV. The output from FILLET is written 
in binary on TAPE 7, which is read by GRADE I! along with the input of 
Table 111. 

The output fro® GRADE 11 is shown for the first case (evaporator ele 
vation of 2. cm) in Table V. Mrst the input parameters and calculated 
parameters are listed. Host of the calculated parameters *«** dearly 
explained by their name. I.e., wick area, etc. We will coment on those 
that require elaboration: 

- 22 - 
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table x - oomrr reow x 


-JSIWHe. F ?LTr**FTAl UAB-WICK HEAT HIP 6 
FILE NA*c VP 9/2c/76 


INPUT VARIABLES A*D FLUID PROPERTIES: 


LICI'TO MU*RFR 

TEHRCPATUPC, 

LI0IIT9 p® MC TTV 

VA P P R DENSITY...’. 

SURFACE TCh-cjrjw 

.lICUIO„VJSPnf TTY 

V A P n P VTSC P S TTY, ...... ........ 

LATENT HFAT. , ....... 

VAFCP PPF^ciipr,,,, 

rHfP^AL "CPMOUCTIVITY CF LIQ... 

SPEC! FTC hcat 

MOLECULAR WFTCMT,, 

FRE r TTVGTrMee PATU p t ,,,,, 

. .61 R AVI IT A T T n w A t ACCELERATION.... 
ht-AT-PJPF GF1 M F TRY, « 


HEAT-PTPF I W e t o t DIAMETER 

WIC* TUTfKHF**, ., 

WICK ........ 

KICK HEIGHT 

WICK rro«-P 0 T A k j TtP 

SPEC I FIFf) ‘■ATI'PATIGN FPACTITN. 
MI N T*”* ALl n l»’A?LE FCROSI1 Y . • • . 

NO. ^eCTFIFn-PORCSITY t>TS.... 
pnenvjry oniNjT nq. ; 

*'I < 5TA» ,,, F. Tu PC IM 

prpn«rTTY A7 PCI NT... 

PCRHSTTy f.c, 2 

DISTANCE TO POINT........ 

PnPO«TTY A T POINT 


110 

• 

2 


fKfclVh 

• 

3 .duCOSE *C 2 

DK-Pft'S KELVIN 

phol 

9 

6.0GA09E +C2 

KC./CU. N 

PhCV 

9 

f» .26329 F+OC 

KG/CU. » 

ST 

m 

1.95290E-C2 

NEWTONS/M 

VISL 

9 

1 .301 30E-GA 

NEW TUN-SEC/SO. M 

VISV 

9 

9.992fc6E-C6 

NEWT lN-SEC/SO. * 

PFC 

* 

1 .1 6^t6E+C6 

JCUL&S/KG 

PV 

* 

1.06C96L+C6 

N/SL. f 

XKL 

a 

5.G932SE-C1 

WATTS/M K 

SHRV 

a 

I.310C01+00 


xnw 

a 

1 .70320E+01 


TP 

a 

1.9i»AAAt4C 2 

DEGREES KfLVIN 

C-ZE 

a 

l.GOCCal 4C0 

STANDARD GRAVITIES 

IGfcOh 

a 

X 


3-GENERAL) 


HP I D 

a 

1 .OOOCOE ♦.CO 

Cf 

WKTH 

a 

A.OUwOCt-Cl 

CM 


AW ■ 
HW » 
OIAF » 

s « 

EPSMIN « 

NEPS • 

MPS * 
• ( S X » 

>EPF ■ 
USX - 


3.B9u61L-Cl 
9.5^25-Tt-Cl 
2 «OC J.w c**C 3 

0. 

0. 


so. CM 

CM 

CM 


J. CM 

r • I •( 1 

3 .OwC‘C>i>t ♦ Cl Cm 

c ..Mai CO f -t 1 
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CAPACITY (WATTS) 



SIMPLE FELT- METAL sLa*-WICK HEAT PIPE 
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WITHOUT EXCESS LIQUID CONTJUUTION 




EVAPORATOR ELEVATION <CM> 


Figure 5. Capacity Vs. Elevation for Metal -Felt Slab-Hick Heat Pipe 
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HW - The wick height Is the height Ir. the heat-pipe cross section 
to the top of the wick from the lowest point in the vapor 
spaces. 

DIAVS - The vapor-space diameter is the hydraulic diameter, four 
times the area divided by the wetted perimeter. 

HVS - The height to the top of the lowest vapor space from the 
lowest point in the vapor spaces. 

PERIM - Active perimeter of vapor space is the perimeter over which 
evaporation or condensation occurs. 

TH1, TH2 - Angular groove feed locations measured positively counter- 
clockwise from the top of the heat pipe [see Reference (2)}. 

FQGRV - Radial input fraction is the fraction of heat input to a 
single vapor space. 

Under the heading "FINAL SOLUTION," the pertinent data on the heat 
pipe at its maximum capacity are printed. In this case, the maximum 
capacity is 31 watts, the fluid charge is 5.7 grams, the vapor Reynolds 
number based on the average vapor velocity u in tne adiabatic section 
and the hydraulic diameter of the vapor space is 260, and the vapor velo 
city head 1/2 c v u^ is 4.8 x 10“ 4 cm of liquid. The radial Reynolds num- 
bers, which are based on the normal velocity of the vapor, negative .men 
towards the condensing surface and positive when away from the evaporat- 
ing surface, is -2.22 in the condenser and +2.22 for the evaporator. 

Next the distribution of stress, static head, porosity, saturation frac- 
tion and static pressure of the vapor are listed. The static head is 
the contribution to stress from the gravitational acceleration. 

The results of several runs for the capacity of this heat pipe as 
a function c* elevation are shown in Figure 5. Two curves both with 
(IFITS * I) and without (I FITS ■ 0) the effect of excess liquid are 
shown. The wick in this heat pipe was intentionally selected to have a 
low permeability in order to show the marked effect excess liquid can 
have on a heat pipe. 

5.2 EXAMPLE OF GRAOEO-POROSITY-WICK DESIGN 

In the second sample problem, we focus on the design of a graded- 
porosity wick. The input is given In Tatle VI and the output is given 


31 
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In Table VII. The heat pipe has an Inside diameter utf 1.1 a# (HPID » 
1.1). a wick thickness of 0.51 cm (HKTH - 0.51) with a fiber diameter of 

0. 0127 on (DIAF » .0127). In this case the wick is horizontal (IGEOM « 
0). The condenser, adiabatic and evaporator lengths are, respectively, 
50, 60 and 30 cm (XQ > 50., 60., 30.). 

The wick is being designed for maxi man capacity when the evaporator 
end is elevated 2 cm (ELEV * 0., 2.). By setting IPRIICD » 0, we are 
requiring the wick to self-prime at the operating tilt. Because fibrous 
wicks exhibit capillary hysteresis, a wick with a higher capacity could 
be designed if the user accepts the operating constraint of first level- 
ing the heat pipe with no load to prime the wick before elevating the 
evaporator end. If this option is elected* IPRIMED is set to 1. 

The TRW version of GRADE II has the provision to automatically plot 
the wick volume density profile (the volume density is one minus the 
porosity). This plotting capability, which is activated by settling IPlOT 

1, utilizes plotting routines outside of standard FORTRAM, and thus it is 
not included in the user's manual. The plot of the volixae-density dis- 
tribution is shown in Figure 6. The wick begins with a porosity of .873 
at the condenser end and ends with a porosity of .631 which is above the 
set minimum of .60 (EPSHItl * .60). The output of Table VII shows that 
the maximum heat-transfer rate is 85 watts, and the fluid charge is 31.1 
grams of ammonia. 
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APPENDIX 

DESCRIPTION AND LISTING OF 
GRADE II 
AND 
FILLET 
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Figure 7. Flow Diagram 
Searches for 


i. 










Figure 8. How Diagram for Subroutine INGRTN, Which 
Integrates Along the Pipe and Reports 
Whether There Is a Failure 
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A.1 DESCRIPTION OF GRADE II 

The structure of GRADE II Is given In Figure 7 where the flow chert 
for the main program Is displayed. The main program begins wltn a call 
to INPT , which reads the data, computes parameters and calls PROPS, 
which calculates the fluid properties. INPT then writes all of this 
information. The main program next calls INGRTN, which integrates the 
differential equations from the condenser to the evaporator end with an 
assumed value Q for the heat-transport rate. If the grooves are found 
to dry up, INGRTN reports this by setting IFAIL . ■ 1. Q Is Increased If 
the heat pipe has failed or decreased If It has not, in a binary search 
for the maximum rate. When the change In Q is less than 1/22, the pro- 
gram is assumed to have converged. A call to WRT writes the final solu- 
tion. 

The structure of subroutine INGRTN Is displayed In Figure 8. IFAIL 
is Initialized to zero and then RUNGE, which makes a single Integration 
step DX along the pipe, is called repeatedly. RUNGE relies on subroutine 
DERIV to supply values of the derivatives of the key variables. When the 
/, Integration Is In a region of evaporation, a call to DRY is made to check 

\ whether the grooves dry up. if they have, IFAIL Is set to unity. INGRTN 

returns control to the main program when the evaporator end is reached. 

i 

A brief description of each subroutine is included in the listing, 
where the subroutines appear alphabetically. 
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000061 2** rnNTT**» ,e 

C 

000061 ran «®t 

C V*T MINTS THE f INAL SCLLTICN 
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A44(6)# 

AAM6J/ 

4.b4711CE>4# 

-A # 90671 AE-7 / 

A54 1 b ) * 

Ai5(6W 

5*803347E-1# 

-3.2S7U86-2/ 

464(6)# 

Afcb(fc)/ 

1 #0 764£GL -b# 

-1.3t3]?36-8/ 

A 74(c) # 

A 7 1> 1 f ) / 

7 c 4 0# 

UlF7bAbf-l/ 

At 4(0)# 

AtbCt )/ 

4*7197lb::-6# 

-S.211675E-9/ 

494(61# 

A 9 5 l M / 

•2#32Lirbk‘42r 

-JU-C7AA*, 9t ♦ 1 / 

4104(6)# 

AltbCC)/ 

1 • ?*4C 7 47*1 >f* t 

-3#V3223*-li/ 
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000005 PATa 4 1 ! 1 ( 6) # 4112(6)# 4113(6)# 4114(6)# 4115(61/ 

000005 *_^^5J77 c 7^r3/-5^463947i-5f 2.856240LW# W*817C0t-10#8*14e495E-13/ 


•V c 





c 

/» 

KITPtiGtN (-3ACF<T<-250F) 



t 

000005 

DATA A J 1 ( 7 ) # A21(7) / 



' 

000005 

* 113. 9# 2».016/ 




C00005 

DATA A3 1 ( 7 ) # A32(7)» 

433(7)# 

A 3A t 7 ) # 

435(7)/ 

000005 

l.AO# o.# 

0*# 

0* # 

C*/ 

000005 

PATA AA 1 ( 7) # mA2(7># 

AA3 (71# 

444(7)1 

445(7)/ 

000005 

* 7 . $a«o7A£+X#-2. 3055566-1# 

5.3175996-3# 

-2.34C715E-5# 

0*/ 

000005 

DATA A*5U7>* A5? t 7 ) » 

453(7 )# 

4 54 ( 7)* 

455(7)/ 

000005 

* 3.?l7!73£-»l#«l.A312e9tTZ# 

3 #C64?64fc *0# 

-3.13377701# 

1 .3 76 A A9 e-2 / 

000005 

OATA AH<7># A62 t 7 ) # 

463(7)# 

AtA ( 7 ) » 

At* (7) / 

000005 ♦ 7. 5»9*71(:E+l#-3. 3232326-1# 

2.281A69E-3, 

-7.A78632E-6# 

0*/ 

000005 

DATA A 71 ( 7 ) » A72 ( 7) # 

473(7)# 

A7A ( 7 ) # 

475(7)/ 

000005 

* ?.102A02L4l>-7.503727£40> 

9*6 J 3273H-1# 

-A.8611 16fc-2» 

0./ 

000005 

OATA • A%2t7># 

433(7)# 

A f At 7 1 # 

A45(7) / 

000005 

* -l.nP7'“S£4l# 3.560937E-1# 

-3.8CC9A3£-3# 

l*6t7G39t-5# 

*2 • 6C9* 66E-F / 

000005 

OaTA A92 ( 7 > # 

492(7)# 

494(7), 

A 9 5 ( 7 > / 

000005 

*' 1.7tP67oE+A#-1.371991c4A» 

4*10394 5£ + 3# 

-5.A53515642# 

2.716t2A£4l/ 

000005 

OATA 6101(71# A 102 (71# 

4103(7)# 

AlCAt?)# 

4105(7)/ 

000005 

♦ 1.17pO0O£-l#-7.992A2AF-5# 

-1.A0X5166-6# 

0 .# 

0 */ 

C00005 

~ OATA Al’ 1(71# A112(7I# 

4113(7)# 

AUAt7># 

4x15(7)/ 

900005 

* l.ft3A'>'Ue-3.-3.768939e-6» 

-4,37937iC~6# 

1.27C-396 r -iO# 

6*/ 

c 

000005 

— TF rt.ro, 01 6C TO 20 




c 

000006 

T2 » T*T 




COO 00 7 

~~ T3 * T**T 




000011 

T 4 • T7*T2 




000013 

tp . l^PC./T 




000015 

TP? • T**U 




000017 

TP* * r QPtTR 




000021 

TP A * TP 7 *1R2 




CCC023 

AIT. At rrtjt 




000027 

AL T5> »At T + AlT 




000031 

At T’.At T ? ♦ A 1 1 
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000033 «LT*»ALT?*4LT2 

c _ . . . ■•■ - 

c’ Pluto PROPERTIES 

C 

O00035_„ TF - 411(1) 

000037 yww«4?l(l) 

C0004 1 SMRV » A'*l(l)*A3<i(L)*T'*A33ll)*T2*A34(l )*T3*A25(l )*T4 

000053 MFC, ■ 441 (L )*A42 (l l*T*A43 ( l )*T2 + A44 <1 )*T 3*A45(l ) *T4 

000063 PV • FjfB(A51(L>*A52<U'*Tk*A53(ll4TR2*A54(L>*TP3*A65(U*TP4) 

000103 RHOL • A61(l)*A&2(l)*T*A63<l)*T2 + A64U. )*T3*A65 (l )*T4 

000115 PHOV __ ■ F*P(A7Hl)*A72(l)'*TR*A73tU4TP2+A74(L)4Tk3 + A75(t>*TR4) 

000133 ~ VT*4l * ■ e *P ( A61 < l ) ♦ A62 ( l ) *T ♦ A63 M )*T2+Af4(L)*T3*A65(t)*T4) 

000151 VISV ■ Pyp4A91(L)*A92(l)*ALT*A93U)*A(.T2*A9A(l)*ALT3*AV5(L)*AtT4l 

000167 **L • Al01(U*A102(L)4T*A103(L>*T2*A134(t)*T3*A103(L)*T4 

000201 sf • Alll(L)*All2(L)4T + A113(l)*T2*A114(l )*T3*A115( U*T4 

000213 VISV • VISV/4. 16975 04 E6 

C$0215 _ VTSL*VISl/4.1697504E6 

0C0217 >ftpom 

C 

C00220 „ _20 C^ntimuF 

I 000220 ~ ’ RETURN' 

000221 FHO 
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«M?oriiTT»Jt RECfcCECR» 0 H,AC» l*EC) 

' c~ * . ■ , 

c rfcfpe Calculates the parameters for a grcove with a 
- - ■ - c circular meniscus 

c* 

000007 /•’ARAM/ EPS»FPSC»PNO#CFf/G,NI # N! 01 ,HW/i:tOM,wKlh,t-PIC# 

000007 J NC»FC{ 13 )*NFL 2 V/XElEUUO)p 2 ELEV(ltl»ElCVUCU 

000007 ? lLeVBI 10 )»H#FVLe»AW»Ae»XTC T#OX» 02 »COTB»COOT»Fi» 

000007 3 OlAF # A{ 6 » 500 J#XfcClCJ# 20 ( 10 )#i>BS» 6 RVS*H: LF,FCCPV# 

000007 4 NCASEf LASTEPS»SAVEPS( 000 >* 1 F A 1 1 PBO» I M« #NH, t PSS , 

000007 ? NVS#AVS»VrF # ClAVS#W#C 6 PTH,PHI»AFG»»tT # S#IFPlMf t» 

CCC 007 ft NEPS*XEPS< 2 O)#ePSX( 2 t»#NOd*SS#Ht ( S»AAA(<iC)# 0 OM<iOJ* 

CCCC 07 __ _ .7 VELHD»ZOP(lC)»XKOB(lC)#XQ 8 (lw)»IflTS»STRS( 40 )/ 

000007 B IFASS»P£PIM, ROUGH, THA,TH 2 #IEV< 1 U#N£V»U#FFM 1 # 

CC 0007 «» IPLOT»XX 0 ,XXl,EPSMlN,hIGH,LOt«,HDl<fc),HO 2 <fc> 

000007 _ _ CO**MMN /PpROPS/ XMW»SHRV»HFG#PV»RHJL,RH 2 V,VISL#VI S V, XK L, SI , T F 

CCCf 07 /r»VOTA/ ACk# PMP, S P, CP# TP# SB AX.CMAX, PS 1 MR, HC, A IP I,f£PIf , 

C 00007 1 RMiMCU»V,BCT,FlQ,Gf AC » Vf AC# R V MN »R PO 

000007 Cr>M*«fiN /AALFA/ Am»BFlT,Cf M, FRF 

“C * : • 

C * 6 AN$ MENISCUS ATI ACHED TO GROOVE UPS 

C IRFC -2 ** 6 ANS MENISCUS RECEDED PROM TIPS 

- C I»E'C*i *fcANS MENISCUS TCUCHING THE GROOVE 30 T 1 CM 

C 

000007 . irfP.LT.RMIM GO TO ICC 

C " ' . 

C •‘PNf^CUS ATTACHEO TC- CROOVfc TTPS 

C 

300012 <®*T»fW/ 1 uO.J/< 2 .*«> 

0 CCC 10 ro«t» r PPT fl ,-SPSI ASPSI I 

0 CC 024 Ptro.ATANiSPSIrCPSIl 

000033 " ACn«.or*,-PHl-<PSIR/RPO» 

000037 f All XLFAIACCM 

000044 *C»f fV/T 9 t.)/ 2 .)Mm»/i 30 .>/( 2 .*TF>)«R*CP$I)-P<lR*R*R 

CC 0057 ™ am •A c —x T® l — 

CC 0061 * ftw«A.*Af /f'tRIF 

CCC 064 TfST*o*f l.-CFSI) 


.8 

I 
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_ SUt>Pt?t'Tr*'fc KLNCtlY.YP.DZ) 

" C 

c OHWGF TAKES CNE INTEGRATION STEP 9Y THE RUHGE-KLTTA METhCO 

C 

000006 ’ Oft'^STHN Y<5),YP(9)»YPI(5)#YP2(9>.YP2(5)#YP4fSJ*Y0(5l 

000006 "fl 1 

000010 VP1ITT.VPUI 

cooois “i vnm-vmtvpi «i)*cz/ 2 . 

0OCO22 CALL *)PPTv( YPi # YD» 

000026 0* 2 1-1.9 ... . 

000030 “» Yntf i • Y ( I >.»P* ( n-OZ/2. 

C0CO37 CALL nrof v < V*>3,Y0) 

COCOAS _ * 1-1 O 

000069 % vom-vti )+YP3(lJ«0Z 

000093 CALL OERTVCYPA.YC) 

000097 PO A T>1 *9 

000061 " A VtT1"YlTW.«VH'in*2.*Vf2ll)*2.*VP3m*VPAm»*t02/6.) 

000076 CAH OPPfVCYf’^O 

000102 PFT!»»m 

CC0103 ' " «N0 



• I ■<&»***■ | * 


» 


I 


•m„c.c»uinfL.ty.?**?*.i!*'v oomm. 13.26.02. 



1 

I 


000067 T»PC«1 

000070 __ T« , n >r ST.*E.(OEHH/lCC.n HP SC-3 

CC0Q75 " teTUPN 

c 

009076 _ 100 . CONTINUE 

C -eNT'Ct 1 ' RfeCE0:0 INTO GRCCVE 

GC0076 IPCC-* 

C 00077 ' »f«w«Avr.vfrl 

000101 : *CAU ALMCACGNJ 

000106 P e ®T»?»*P*SfAX/SP-BCT 

0001)1 AC»P*S«A 7 *(B*SHAX/ 1 P»P*C»'AX»-P 5 IPR*R*R 

000123 ' Af«*f-ATB1 

S00126 OM.t.*Af /PERI 


000130 TF«T1 .0** max/TP-HG 

000)34 ' r?str«p*ei*-cfiAx» 

000117 IMTF5T2.GE.TEST1) IREC-3 

CCC1A3 P^TI/pn 

0001 A A r-o 
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000002 
800002 
000002 
C0C002 
000002 
C000C2 
OS 0002 
000002 
000002 
900002 
000002 
000002 
000002 
000092 
0*0002 
000002 


•Jfopni'TTMfc VS9* 5 _ ' - 

V'AKR CAtCUlAttS f Ht VABCR-SPACC OAR $m$s» WHICH sirs 
TMF INITIAL STPtSS POO 

rr*rr» /® ARAH/ t^S#fPS0#^0,C£f #6i.Nl#^lPl#H,,ICeOH,.RlH#>J 1C, 

I NC,fCaCJ,^tl£V»*6l='-C10lf2ktCVtl0»#ftlVllC I* 

• i *. fVM10t»H,f M'P, AW. AB.HTLT »D*» 0 2 *001 o# 0001 »H)» 

I OIAf, A(f#!.CO>,>OUCJ,2C»10>,*»8!#6PVS,HAEF#f CCt-V, 

«, . ScA$t,tA5TtOS#$AVr0SI900»,IfAll#FB0,IHR,hH,EfSS* 

— . KVifAVS.VJffOlAVS.W.OfcRlH.PHWARtWkt.S.tBAl^kt* 

1 Hf PS.Xf *SUOI#iPS*i2t.».NJP#SS,HVS,AAA«i.U#CCM9v)» 

7 WELHr»2CBU0),*KO8UO)»>«Jt»U0»#lHTS,S1RSl‘O># 

-- m IAASS»PfHM»RCtOH»THl»TH2#IfcV(10J»NeV#PI»l H*l# 

e IPL01,AX0#> «l»?PSPIN#HlOrt,lO»#Hf 1I6J »H02C0> 

CO«f«PM /C PROPS/ XHW,SHBV,Mf6.PV#PH0t ,RH0V#VI5L,VISV#XKl#ST#l 

ojnew^IPK $ a ( 3 ) * SB 1 3 I » SO 1 31 » SM 3 1 
QAtA ( S A 1 1 ) # 1 • 1 » 3 > /.6«» J.‘j#2.2997/# 

1 (3* Ml, 1*1*3) /* 9 C».aO» *S0/» 

— j»‘ " t$r(I)»I*l,3) 

3 f*Hf 1 )#]•*, 3) /i.9#1.2#l./ 

T ? (crf.eo.O. I CO 1C 10 

• •P.<«T/<HPIL/2C0»>»IRHCL-PHCVIA0£c*0*Hi«/190.»/PNj 

«»1 , /I 9 « *St / I (RHOC-R hCV >*G fc£ *G* €01 AVS/190. 1**2 l * 
TTMi?/»9«SA<J»/eP*S0CJ» IMl.-$NUI*6*P<-S6C Jl'SCRMtni 

YT«TT9tf»lAVS/lC0. „ „ 

- **ftt.(YT*(Mb-HVSI/10C.J*< e HCl-XHCV)*C-.C*C/?H0 

im»«m t » peoi ) p»w«ppci 

Brtl'PV 

10 #«P*14,*SI/I0IAVS/100.11/P**0 
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WM COHHlE» 


c 

c 

c 

800010 
OCCCIO . . 

QC0010 
000010 
000010 
800010 
QCC01Q 
CC0010 
000010 
00001 0 
OCOOIO 
00C010 
0CC010 
000010 

' ~'c 

800010 
000013 __ 
ecoois 

COCCI 7 
0CGt/2c 
000025 
CC0027 
CCC031 
CCC039 
000091 
8CC095 
CC0OA7 
COCO 51 
C00053 
C0CC5 V 
(.11057 
OOCOt 1 


< 


«t»90r*iTIv t h* TCT$TCN»R1»R2* U1.TL2I 

VfT fHT^ORATES AICKC 1HF ORCCVt' FROM THl ANC TH2 ANO 
Of OORTS IOC Allows Til ANC U2 ««€«£ ORY-OR OCCURS 


t 

2 

3 

A 

5 

6 
7 
* 
9 


1 


(•OMVPM /r*R4 1>/ t'-Sfl <*SC»RNO»GFI *6»’»l»NlPl»K**ICtOK#WKTM»HPIO* 
NS,rC<tC)**CUV*X£l5VllC)»2Et £VliO)*£lEVl ID* 
£LtVe(30)»H»FNtB>*WtA8» XT CT #D> * C2 »CDT8» CO Cl * R6» 
O;AF,A(e.»OO)*XC(lOI»Z3(lO>*B0S*GRVS»HRtF*f CCRV* 

NC ASE#l ASTERS* SAVERS* 5001 * l FAIL# 280# IHA#NM*fRSS* 

*OlAVS#t**t£RlH*RHl*AK».*T*S»lPwll'tC# 

Net'S* *tRM2C)*f RSM20)»N0?» ST#hvS*AAAI<.U) ,LCM*0)* 
VcLHCi/rRIlOliAfCBI lUI,*0M10».irtTS,STRsl9CI# 
IRaSS* b FR1^/'’CUCH*Th 1*TH2» i£VllGI»NEV»PI»fTI'l* 
IRUOT* AkO»x*l *EfSKlN* HIGH* LCWfKt 1(6) *HC2lfe) 
frr*r»»| /rPRCi'S/ X«te#SHRV#Hf 6#PV*RM3l*RH0V#VISt#VlSV#XKt» 5T*TF 
cm*««rv /r.rtvct// ACR.FM* »SR»CI»*Tf »S> a*»CHAX.P$ IRR.WG* AIR 2 *HPI( * 
HlN*CuRV* 0 CT*flC»'»f AC»wrAC#RWrN*RRD 
f /*AUA/ AFIT*0F1 T*CiIT»IRF 


NTJ»tST««-THl 
TRINTl.eO.OI NT1*1 
*Tl* w Tt 

nri»(T<Trh-fH)/«i • 

«TR«TM?-TSTCN 

TFINT?,fO,0)NT2«i 

*T?***T? 

OT?«(TM?-( ittNJ/ita 

isic(w-SHil/ITH2-TN2 I 
FI n* < r*2-19fOW)/( Im 2-TN1 ) 
rF('*i«rt r, ;/xU 
nrt«v.rt»'*/»T* 

•1*1. /•! 

e.n 

rri 
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3CCC63 

CGCC65 PT.ptl ’ . ■ 

000067 7»T**1 

000071 m*6»0 

000072 100 CCNTtHU' 

000072’ '*<' 101 !•!#* 

€00074 *•!./•* 

000076 C*ll SfffOf <*t#CH»4C#l*fC» 

000104 ” irf f*rr,fc.ll GC TC 102 

0001C6 0*rT.-cr*C*SlNC6P0*n*fPf*Vf AC*FF10/UC*0M**2I 

09012 1 0®ni«O»OT*«PO. 

000127 

000132 reir.cri «-of tc 

000134 T.T46T 

000136 101 cnNTTNlir 

C 

000141 102 CHNTTVUt 

000141 ‘ 1*1 *C»!U AOM 

000143 1* f f f t * r .46.2 1 CC 1C 103 

000145 m-t 

000146 H-4T2 

C001 50 «*•»? 

000192 FHP.-nPi' 

000154* ’ pFl»'«-f‘riC2 

000156 f*T«»0t7 

0001 60 T*TM? 

000162 *■* 0" T« 104 

900163 lOf CPWTINl'f 

000161 tl?-T 

000164 * **1l'»*» 

000165 *»«* 



*-*«**-• **gr» . *«* SjmF -t' 


-V 


RtNX .(Off 111 * jvr « ,?,?«•» 


C4/2F-/76. 13. 2%. 02. 


r.i- 


j - 

*-i 


|: • I 

% * 


- i 


300002 
300002 
300002 
300012 
000002 
300002 
300002 
30CC02 
J0C007 
300002 
300002 
300002 
C 000 02 
300003 
000014 
00001* 
900024 
000031 
90 00 *0 
9CC043 
0600*9 
0(0047 
900091 
000064 
CC0067 
t 6(071 
CCC101 
OOOill 
C00114 
6(0121 
CCC123 
CC01 ? 3 
000121 


c 

c 

c 


f *rt hfti 

OPT WP|Tes TNt atSUlTS 

• * ?«**:>»! /Pit * nt E?*»£PSG»f *0#C( f ,6*41 .NIP 1»H«» Icf Cr»-Kth, t-r IP, 

1 ».C,fC<10».»iSU WHflfVllw! WlUVI101#fUVtlC># 

P ELtV8(10l,h»f fit e,A*»Als**UI*0*#C2#3GTto#30C;TjMs» 

9 0IAF,A(P,*,C,))»»C1 10»#2C,I10>#PIS*C»VS»HR£F # »CCPV» 

4 NCASF.l ASIEPSp 4*VEPSI5031#ir*Il,PR0,lMh,NH,tHS» 

9 HVJ .AVS.VM »CUVS#*.»rtPTH,Pril *4NCbtI/S» IP* IK 0# 

4 .it* $ v »tfrSI2wf»€FSxt26lrHJP*SStHlS**A*t6w)-«CCf(ftC>» 

y Vllh3t2r,P il ;i,»p 081161 »*CS Oil » IF HS.Sieil 401# 

• 1 pals, fit- ir#PC10H»m»lH2, UVUC),Ntv,K»M f-l, 

"o lPtCT,X*0»>*l#EFSPn,hlOH,L0rf»M0H6l .M02(fc) 

coupon /'“ Props # *« i .# SKPv » Mic # Pv # PMni *« HO %, viH , visv #» Pt » si # i » 

n » l * r »»< jn * Kt Vt ( 10 ) * a ( 11 ), w«(l 1 1 , »H 11 1 

. r***6^**«(6#l»!7.1 IMAw/luCU .)*P6U*(*ttT/lCO.I*UOO. 

<*hT»n«T»*tOC1 

VVAP«(C0T/HFt)/(fL0AT(XVS)«(AVS*l.£-4)*RH0V) 

prv*VVAP*(01AVS/100.l/(VlSV/PFOV) 

VflMO»100.ft(.l*RhOV*VVAF**21/((RHCL-SHOVl*6» 

TP fPP» .l T.PCF H I CL (C ti 
*(1l«4. 

(*" **1 T • 1 » 1 w 

* » W ( f » 

p7»?.*Aip-cU((2.9i*p(ii/ptv«i*fltrH/n. ?«ot avsi » 

WW(T1«»1**2 
|r«T..P) /*2 

VP(»|.l.*(>.«2/PtVI/(2.11‘*(l)/KlVP8Ll.t»H/t3. ?K JAVSII 
*»(♦?»•»( l-fchUI/kPUl 
F»|» /*((♦!)••£ 

TF(AP«ITf Sl-i.l .11. .(Ul CC to tl 

91 . 

4 i r** M i **• r 

0 * “P l * 1 #XK ' 

V«’*rf . i r c 1 1 1 •wfl /<i» f |/(M«CV*( *llll/lvP.)*(4Vi*l iM-/;C3. 11 


. 

30 

*r 
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09 / 28 / 76 . 13 . 29 . 02 *. 


800190 *?V»# II*VVAPR.«0XA*.$/2C( • )/l V HV/BHlV) 

8001*7 79 

000152 ;T ““ OH *0 J-l.Mfl 

000199 AM. J»*»T0M*<1M> 

000160 

000163 TrMn.jt.ct.mEviiM'n i***im*i 

000172 M* .31 MHt VI. M*!)* (Ml* 3 l-mCV(lPK”lll»tcUWlM«)-tUVtlrK-ll) 

000172 1 /<*6UV(I*x»«xHtVf IP«<»21)M6tt 

900210 trtrrc ,* flt (,.) *<6.31*0. 

•00211 SAVMIJI* AI»mI 

OOOtlO ... AI«.it*AI«,3l*Mi 

000221 9* 

000229 W»!tr Ifc.CfcOI 

000230 v»!t« t*,*>t»21 Ctl#iN*SS»s t V,7t t K> 

600298 I*f**V .01. 23CC.) 861T((6»967I f 

000299 K*ITrf6»«6ll 

000260 oe *0 9*t •«•€ 

800262 """" v»|Tr#*,Ofej| l,«fy*fll 

088272 60 fn**?IN«>.- 

000279 v*ftr t**«»6*| 

aiooof - •>" 9f> /•i,Nifi 

OOOMt lit 1 7716.9661 A(1.3).A(2»3t»A(6.31.9(5.3l»A(?»3l»A(9.3t 

Off M2 

000119 960 MMMMui. + MlSMl 1*91 SttlllC* //I 

000399 96* rpe«.Tn0a.4w>*lHk P9*I*l.P HCAl-TMuSJ t* «9Te IS... ... . 

000319 1 *1«.9»7M 6*115/ 

000319 — 10*,9ChlHC 1CTA1 MCUtO IN *ICK IS. ........... . 

000339 * *12.9.71. «**«S/ 

000119 8 10>.9vHlhC V 9*0* 7179810$ NurttA IS........... . 

OOC339 9 *12.*./ 

000339 6 1'*»,9v*«1n( HAS. VMC» VtlOCtlV h£90 tS. ....... . 

86033! 7 r 1 2.9.9I- CP UC.l 

000139 " 96» *P6f**r(1(ti # *6HlHE PA*. 11**861967 MICTION MCK» IS*. . 

060939 1 *12.9/1 

000999 9*1 r*»*MTM'*9#13HlHC 9A0M1 *£7KC10S NIMirfS A# (I 1 

60039* * ”“961 WMI1**.ltHiW Ss.CU!* NC. .12.17*. » 

000339 1 *!<.*.} 

8CC339 9f 9 ***** ti/M 




03-^8-9209 -C9292 








09m/76« I 


1 |»j # 'KHlSUikCCt £54 #Ti»11HS?AI 1C 

t ' W(C UC«I»61>9H(CII 1 10« I#36«#9HICH 110*I//I 

.•ITU** ... - > • - 

*mc 






Mm •wM-tmi 
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A. 2 PROGRA* f’ltl.CT 

The mathematical methods employed by FliUT ire described in Refer- 
ence (2). The first pert of the progress computes * separate table, for 
each type of fillet or puddle that can enist in the heat pipe, of Jx 
area, free pe rice ter, wetted perimeter and hydraulic diawiter for various 
values of stress. By statement 500 this is accomplished. The rest of 
the program manipulates the data into a usable for®. For the fillet 
that form at the bottom of a vertical wick and the puddle, the values 
of stress in the tables does not increase aonotonically. SiA>ruutine 
RFARNG rearr&ngev these tables for iiscreasing stress. 

The next step is to obtain total values for the area and hydraulic 
diameter for all fillets and puddles in the heat pipe. For a specified 
value of str*;s, subroutine t.%T£R Interpolates the tables, and the a«jw$ 
and wetted perimeters are sirred froa which a total hydraulic diameter 
is calculated. This is done for a range of values of stress to construct 
a table of total are* and hydraulic diameter as a tunctien of stress.) 
This table is then written on TAPE 7. 
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09 / 28 / 76 . 13 . 05 , 23 . 


June* A*. F ?UtT ( IN f U T ,OL T F 0 7 » T AP E 5^2 N P 0 T> TA P fc t » OUT PUT / T A F £ 7 ) 

PUIFT CALCULATES THE MENISCUS SHAPE# ARIA, FREE PERIMETER, 
WETTED PERIMETER AND HYORAULIC DIAMETER. FQR_. FILLETS IN A _ 
SLAP -WICK' HEAT PIPE AND A BILGE IN A CIRCULAR TUBE. 

NAMFLIST /f ILL E TO/ HP ID, WK TH, RHO, ST, 6c £, IGEOM, T KE L VN, 1 10 

R IS THE TUBE RADIUS IN CM 

T IS_ TH P HICK THICKNESS IN CH . 

RHO * IS"THI DENSITY DIFFERENCE BETWEEN LIQUID ANO VAPOR 
IN «G/M**3 

ST IS THE SURFACE TENSICN IN N/K 

GEF~T$ THE GRAVITATIONAL ACCELERATION IN STANDARD GRAVITIES 
LIO IS THE LIQUID PARAMETER 
TKFLVN IS THE TEHPERATLR E 


ICNFG 


FOR A VERTICAL WICK 
FOP A HORIZONTAL WICK 
FGR NO WICK (BILGE ONLY) 


DIMENSION $TRESS1(80),STR£SS2(40>,STRESS3( AO),STRESSMAO>, . 

1 S1RE$$5<A0),STFE SS6(P6),SF1(B0),SF2(AC),SF3(A0), 

2 SFMA0>,SF5(A0),SF6<e0>.01<60>,D2<AC>»D3lA0>»DMAC>, 

3 D5(*G),D6ieG>»AlieG),A2(*0),A3(VO>,AA(AC),A5«,C), 

~A ‘ A6(8C),SWi(3C),SW2(A0),SW3(4O),SWA(AC),SW5(A0),SWfc(80), 

5 DDH(tO),AAA(FO),STRS(80) 

■“enWR’ /f PROPS/ XHw, SHR V ,HFG, P V, RHJl ,RHOV, VI SL , VIS V , XhL, ST, Tf 

PIHPNSIHNLESS parameters 


PR AD FILL ETD 
TF(LTQ.EO.O) GO TO 10 
' T R A R'P » TK F L V N * 1 ,8 
CALI PP n P SC L 10, TP ANK ) 
o,un* .2 31** 3* < PHQL-PHCV ) 
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A 


UNX COMPILE* fVgR.2.3 Hj 


09/26/76. 13.05*23... 


FILLET 


00020 


ST-6. 66**3. 2616ST 


00022 

00026 

00026. 

000*30 

00037 

0GO61 


0006 S 
00066 
00051 


00053 


TfS G-9.P0 

ICNFS-IGPOM 

TF(JGEOM.eQ.O) ICNFG-2 
OPT T / < RHC*GE S * G ) ) 

R-HPIO/2. 

T-W*TH __ 

PB»tR/100. >/A 
TB-!T/100.)/A 

PI -3.161 592656 _ 

OaTpha-fI / lbl. 


C. 

C 


SET RANG E fo r TH|.,S TOE $ S . H „ AT THf INTEGPAT 1CN ..STARTING .POINT. 


00055 

00060 

00062 

00066 


C 

c 

c 


•00075 

.00077 


C 

C 


H0*l./tl0.*RB) 

HF-?Oj[ . . ..... .. „ 

Vh-60. " 

GG- 1 HP /HO) 6*.( 1./ (MH-1) ) 

" CALC tO TIO NO FARE A / «ETtF0'>"ERiHeT6R» FREE PERIMETER/ 

hyopaulic oiameter and STRESS at each value of H. 
on' - * ob~v -l * nh ~ 

H«H06FG6*(K-1) 

TOP"FtLLeT WICK VERTICAL (CASE 3) 


(00106 

foouo 

.00112 

100116_ 

(00116 

>00120 

(00122 

>00126 

(00125 

(00126 


TFttCNFG.NE.il GO TO 66 

y.-f 3/?. 

AL®HA» p I/2. 

Y-H 
XO-X * 

YO-Y 

ALPHAO- At FHA 
AA*0* 

sr«c. 

f>n 


TS© 

►o o 

gg 

- 

c 5 

r 2 
-* S3 


vD 

on 
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I 




jttUU&gBMiJH t . (V 6B f 


09/26/76. ,13.05.23.. 


HLLFT 


CAU STMPSN(X*Y*AA>SF, ALPHA, 0AIPHA.X0.Y0, ALPhAC 1 

CRIT-X-sn*SlM ALPHA J 
" IF ( C Pi T .GE. C.J GC TG 22 

70 CONTINUE 

c-n to_? 4 • ©fi - 

Y? VO— Tn/2. *** 3 

Si 

W 2 *x . .. ... o£' 

W3*wc 50 r 1 

VO-M 

VI «Y*RB*COS (ALPHA) . . . . - - £ > * 

V?.Y 

' V3-Vl + $OPT«R6*RB-wO*wC)- v-3 ,- 

t-i, CO 

HfHINAT’E UNSTABLE SOLUTIONS 

If (STPESS3CK— 1>*N6»G» .AND. V1.GT..STRESS3CK-1)) GO TG 26 
STPP?S3(N-1)«0. 

A3(*-l)»0. 

«F3t*-l)»0. . - - 

YV3(K-1)«0. 

03<K-1)«O» 

26_CALl $ANf CA3{K^?.W,W0*Wl#W2,W3>V0,Vl»V2*V3,AA) 

«TPE C S3(N MV1 

$F*(K)»SF 

$V? (K ) .PW+V3-VC 

n?»V)«A.*A3CK)/SW3(K) 

tnwFP AND UPPER BOTTOM FILLETS — MCK VERTICAL (CASti 1 * 2> 

?A Y-TO/?. 

ALPHA— *>T/2. 

Y«H 
XO •* 

YO*T 

IF (Yf'.LF .1.0' GO TO AA 
AL ,, VA'*»At CHA 
A A • 
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UNX CO H P 1L6R ( VR R.2.3H) 


09/28/76. 13.05.23. 


SF«0. 

TfLAC-0 _ ■ ■ 

00" ?0 1*1# 360 ... 

CALt S IMP SN < X# T»AA#SF» ALPHA# DALPHA#XO» TO# ALPhAO) 

CR1T»X-R8*SIN( ALPHA). . 

YFUFTAG.EQ.1) 60 TO 2t ■ . 

ITICRIT .LE. 0. .AND. ALPHA .LT. PI/2.) 60 TO 

IMCPTT .GT. 0. .AND. 1 .FO. ItO) GO TO AA 

TF fCRIT~« GT . 6. .AND. ALPHA »GE. PI/2.) GO TO 2A 


TFICRIT .GT. 0. .AND. ALPHA .Gc. Fife*) bu i 
CONTINUE 

CONTINUE . . . . - - 

WO-TR/2. 

Wl *0. 

V?»X - • 

W3»W0 

VO-H 

Vl«t*PB*COS(ALPHA) 

V?*Y 

V3»V1-S0PT (RB+RB-WO*WO) 

CALL SANF ( AA1#PW# W0# hl»W2# W3# VC# VI# V2» V3# AA) 
kl»K 

IFUFLAG .60. 1) Kl-K+NH 


IFUFLAG .60. 1) Kl-K+NH 
IP( IFLAG.EQ.l .AND. 101. E 
IF { I FLAG .6 6.1 .AND. STRES 
TF( IFLAG.EO.l .AND. VI. GT 
TrtiriAG.EQ.l .AND. IQl.E 

CONTINIJF 

A1CKD-AA1 

STPFFSKKl )«V1 
*V1 (K1 )*<> W ..VC-V3 

«F1(Y1 

P1(X1)*A,*A1(K1)/S’«1(K1) 

T F C I F L A R .60. 1) SO TO AA 
iflap*i 

GO TO ?<? 


I01.E0.1) GO TO AA 
STRESS l (Kl-l)»fU. f i») St. TC 35 
Vl.GT.STRESSl(Kl-l)) 101-1 
IQl.E 0.1) GO TO AA . 


;rp FT|.Li.T — HORIZCNAL WICK (CASE A) 


vo 

i 


/ . 

FILLET — - - — 




t — 


UNX COH RU.ES-. lyaUMfli- — — - 09/2P/76... 13...05.23., - 

00510 44 IFMCMPC-.NI.Z) 60 TO 68 

~ — ' ' 

.00513 7«H 

>00515 ALPHA*PI 

tmin^ - " 

iAOH?! Y0«Y 

IC0523 1HY0«IF «E* ) 60 TO 56 

100526 - alphao-alph* - . 

100530 * " " ~ AA«0* 

>00531 SF«0. 

^C0534 CttCL^^T^PSN? X* Y#AA# SF p AL PM A»0 ALPHA# X0# Y0» ALPHAO ) 

IoSsaJ C R TT.R^C0S(AUHA>4Y-H-rB/2. 

300555 IFCCOIT •*6.6. 9*> 60 TO 52_, 

500557 ** 30 CONTINUE 

000561 60 TO 5A 

000562 32 W0«0*. , • • - - ' * ' 

000563 .' Wl-»-*3*3 IM ALPHA ) 

9003TJ ; «!• • 1 • — -• • 

000601 V0*H 

000603 Vl»H+T8/2* 

000606 Y?*Y. . - • * * 

JJJJK GAIL SAVF(A4<K>,PW,W0,Wl,W2,W3,V0,Vl,V2,V3,AA> 

000625 ' ; .\JJ*****<*l m M -- - - 

000627 ’ SF4(K>«SF 

000631 * SWACK>»PV**WC-to3> 

000635 04 (K > *4**A4 (K ) /Sw4 (K ) 

£ t'popp rjLlET •” WICK HORIZONTAL (CASE 5) 

C .... 

000641 * 3P X"0. 

000642 7«H 

OC0644 ALPHA *0» 

C 00645 " V °“ v ‘ - 

000647 

CC0651 ALPHAN.AI ?HA 


s 

VO 

an * 











RU NX COMPILER CVE R.2.3 *) ... .09/28/76. 13.05.23. ..... 

000653 4A-0. 

000654 _ ... - - ■■ ' — • 

000657 CALL^SIHPSNtX# Y#AA»SF» ALPHA# 0 ALPHA# XO# YO# ALPHAO) 

000670 c m ^V+WCOS C ALPHA 1-H4TB/2 . 

000701 * IFtCRTT »LE. 0.) GO TO 62 

000703 60 CnMTINU B 

000705 GOJTT! 6 8 ... . ... — 

000706 6? WO-O. 

000707 W1*X— R8*S IHC ALPHA) 

000714 • W2»X . 

000716" W>~-Wl«-S0PT(RB*PB-TB*TB/4. ) 

000725 VO-M 

000727 Vl«H»TP/2. . ........ - - 

000712 ' V2-Y 

000736 CALL SANF(A5<K)#PW#W0#Wl#W2#W3»V0#Vl#V2#V3#AA) ..... 

C 007 51 ~ " STRF«55(K)«V1 

000753 SF5(*>-SF 

000755 <tW5(K_)«°W + W3-W0 . . 

600761 DM * ) *4. * A 5 ( K)/ SW 5 ( K ) 

C 

C BItGF — NO WICK. (CASE. 6 A .7) . . 

c - 

000765 6* X«0. 

000766 ALRHA«f>._ 

600767"'""" Y«H 

000771 XO-X 

000773 _ 70«Y 

000775 " ALPHAO* ALPHA 

000777 AA*0. 

ooiooo s F »n. 

001002 IFU?m.*YC> .IE. PB .AND. 1FLAG .tO. 0) 1FL AG*1 

6CilOi6 nn TO 7*l>16o 

001020 " CALL ST** t *SN(X‘»Y#AA»$F#ALPHA»OALPHA»XG»YO»ALPHAC) 

001031 CRI T *X-R l> *ilN(ALPHA) 

001036 ..• IF(TFLAG.tC.l) GG TG 66 


LO 




MINX 


09/ 26/7 6 ♦ 13. 05 .23. 


ALPHA .LT. PI/2.) GC TO 79 
1 .10. <JC) GO TO 89 
ALPHA »6T . PI/2.) GO TO 79 


IPCCPTT »LE« C. .AND. ALPHA .LT. PI/2.) 6C T 

IFtCPTT .Oil 0. .AND. ALPHA .GT. PI/2.) GO T 
CONTINUE 

CONTINUE — 

~wo*o. 

wl»0. 

W2- X - 

W3<*0. 

VO»H 

V1«Y*RB »C0S ( ALPHA)- ... - - 

~V’«Y~ 

C A L l ^ ^ A MF (AAG i ^ P- - ®.(L?l ijf if. ^ W».V 2» V3» AA) 
K1*K ” 

TPCTFLA6.60.1) K1«K+NH 

IF ( IFLAG. £Q .1 .AND. . I06.E0.1) GO TO 89 . 

iFciruGiEci.i .and. s mss6 k w).eo.o. go 

IFUFLAG.EQ.l .ANO. VI .GT ♦ >TRE S$6(K1-1 ) ) IC6 

IF ( IFL AG^6Q*I.. *ANO. IC6.E0.il .GO 7^ 6<( 

r CONTINUE 

A6|K1)«2.*AA6 

ST®ESS6(_Kl)*yi. 

?W6(K) )«2.*PW 
?r6<Ki)«2.*SP 

06(K1)«9.*_A6(K1)/SW6CK1). 

— TF (IFLAG. EQ.l ) GO TO 69 c 

A IFLAC-»l •». 

60 TO Hs 

9' CONTTHUP § 

n continue tc 

CALL RPA»NG(STPESS1.NH) 0 

CALL PEAPNG(SF1»NH) c , 

CALL PFAPNG(SVl.NH) £ < 

CALI 0FA8NG( Al#NH) ^ 

CALL PFAPNGtSTPES jfc»NH) kJ ► 

CALL PrA»NG(SFfc»NH) 

CALL Pf A°nG( S1k6#NH) 


GO TO 
IC6-1 


o 


ORIGINAL 

Art — . 


/ 

FILLET. . - 






f 


/ 


A. 


RUN X C OMPI LER < V ER«g»3N) 09/26/76,. 13,05,23, . .FILIJET 

001246 CALL RPARNG(A6,NH> 

001250 n* ?*NH 

001252" DH 600"I«1#NH 

001254 . STRS(I}-H0*GG*4(I-1> 

001263 IF(ICNFG,NE.l) GO TO 507 

001265 CALL INTPR( A3» STRESS 3* NHfSTRSI X)»YA3) 

001271 CALI INTPR(SW3»STRESS3»NH,STRSCn>YSV.3) 

0012 75 CALL I»TER{Al#STR£$$l#N>$TRSm»YAl) 

00i301 1 CALL INTER (Shl;STRESSl*N»STRS(I)»YSWl) 

301305 *. AAAf I>-?.*(YA1+YA3)*A*A 

001312 IP ( AAA( I ) «EQ,0« ) GO TO 600 

001314 "pSW-?.*TySw3+YSW1)*'a ' 

001320 OOH(I)«A,*AAAC I) /SSW 

0JJ1324 _ STR5(I)«5TRS( I J*A 

001327" " "’ GO TO 600" 

001330 907 IF ( ICNFG « E0.3 ) GO TO 517 

001332 _ CALL INTPR<A4,$TRESS4#NH>$TR$m#YA4) 

001336 * "CALL INT e R ( SW4* SIRE SS4 j NH* STRS ( I ) » YS h4 )" 

001342 CALL TMTFR<A5»STRESS5>NH,$TR$m#Y45) 

0C13 46 CALL TNTPR<SW5#STReSS5»NH,STR$<n#YSw5> _ 

001352 917 CALL INTPR(A6* STRESS6»N*STRS<n» YA6> 

001356 CALL TNTER(SW6>STRESS6«N#$TRS(1T#YSW6) 

001362 ■ AA A f jf |» ( 2, * < YA4 + YA5 )♦ YA6 ) * A*A 

001370 '" * f p ( AAA ( I ) , EC«G « ) GO TO 600 

001372 $SW»f ?.*( YSW4 + Y$k*5)4YSWt )*A 

ooi3?7 pr»m t)-4.*aaa( d/ssu 

001403 STRSm-STRSd)** 

001406 600 CONTINUE 

_c 

001411 WRTTC«7) NH*(STRS(n#I«l»NH)#UAA(!) # I»l»NH)',<COM( 1)',I-1,AH) 

C 

_ C TMC Ol'TRIIT IS A TABLE OT AREA IN K*42 AND HYOP AL'LIC 

" C DIAMETER IN M AS A FUNCTION OF STRESS AT THE TLBE CENTER 

C IN N OF LIQUID. 

C 

061441"'"" ” “ ‘ STO»* 

001443 PNO 


I 

o 
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A 


/ 


\ 


/- 


ftUH?L-CO H . p IJ-.EP-A ¥f 


09 /2E/76* 13.05 #2.3* 


0000 1 <> 
0000^0 

000050 

COCO 5 2 
CC0075 
0C0O76_„. 


0X*CnSULPHA>/<2.*7) 

RETUPN - - 

FNO „, „ - • - ' 


© o 

-•*3 » 

►o cS 

gg 

c g 

C=3 

^53 



O ... ,, 
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ftUNX CO V ER. 2.3* ) 


09 / 28 / 76 . . 13 . 05 . 23 _, ... 


900010 


c 

c 

.c 


. I NT 68 l F, X* N# XX, Y ) ..... - 

SUBROUTINE INTER INTERPOLATES TH£ CALCULAT IONFO OATA 
oi^pWroN"F*i),xa» 

XMTN-I.F 6 .... .- ... - — - - - 

'VMAX«-i .E 6 
DO AO I«t#N 

JFixm.FO.O.) GO . TO. AO 

IfTvU ) .Gfc.XNIM GO TC 30 
yt»I v»v ( I ) 

TMTN-T . .. . . - - .... 

' IFixrn.LE.XMAXJ GO TC AO Q 

XNAX-XCI) SfW 

IHAX-I . hjgS 

"crwTiNUF og 

TF(YX.GT.XMIN) GO TO 50 gj j> 

T-.FUNTN) - jn ** 

if TlltH <— i ’’O 

IMXX.LT.XMAX) GO TO 60 o 

.Y- 0 . 5 W 

RFTUPM h< 53 

CONTINUE 

_/?•<*♦ 

X?i«. 

on 70 i« 1 ,n 

IPfxm.FQ.O.) GO TO 70 
F!«r? 

r?«F(TJ » 

Xl-X? 

X? •XIII 

TFCX2.GT.XX .AND. XI. IT. XX) GO TO 00 
CrtNTINNF 

r Y«Fl4.IF**-ri)*T>X-Xl)/IX2-Xl) 

P C TIIRN 
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RUNX COMPILER tVfa» 2 » 3 N) . 09 / 26 / 76 . X 3 .U 5 . 23 . 


_ni«»rUTI»'E RE APNG( A»NH) 

" C" 

C RFARNP rearranges the parameters for the bilge ano 

_ c „ BOTTOM FILLET/VERTICAL WICK FOR INCREASING STRESS 
C 


000005 

DIMENSION A { t-G ) * DUMMY (FO 1 

0COCO5 

no in t«i,nh 

000007 


00001 I 

0UMMY(n«A(2ANH-I + l) 

0000 l A 

OUHWYiI\)«A(U 

CCC0X7 

10 CHVTIVUF 

C00022 

N«NH*P 

000023 

no 90 i«i,n 

000025 

AC j)»bUMNY(I ) 

000030 

20 CONTINUE 

000033 

RETURN 

0C003A 

~ ’ FMn 
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lion. £QHeiLi9..XYFj .urn. 


C 9/2F /76. 13,0- .23. 


000005 

OC0005 

000005 

000005 

000005 

000005. 

000005 

C00005 

000005_ 

000005 

000005 


c 

c 

JL 


1V££9KU!!& f ROf-S U<t 1 1. 

TMIS ROUTINE COMPUTES FLUID PROPERTIES PROP DATA TITS 


COMNHN /CPROPS/ XMW,SHRV,HFG,PV,RH3L,RHPV,VISl,VISV,XKL,ST 
OIMFNSinN All (?) » A2K7), 

A34<7), 

A44{7), 

A54 ( 7 ) , 

A64<7), 

A74 ( 7 ) , 

A64<?)» 

AS4{7>, 


,7F 


mm* 

A41(?f # 
A51C7)# 
_A61(7), 
A71 C 7) * 
API (7), 
A«l(7)» 
AlOi 1 71# 
auk 7 ># 


A32 ( 7 )* 
A42 (7)* 
A52 < 7 ) # 
A62 (7 )» 
A72(7). 
A82 (?)» 
A92<?) » 
A102 ( 7 ) * 
A1 12 ( 7 ) * 


A 33 1 7 ) * 

A 4 3 ( 7 ) * 

A5 3 ( 7 ) # 

A63 ( 7 ) » 

A 7 3 ( ? ) » 

A 8 3 ( 7 ) * 
A<53<7)# 
A103 ( 7 ) # 
Al 13 ( 7 ) , 


A 1 0 4 ( 7 ) » 
A114 ( 7 ) » 


A 3 5 ( 7 ) # 

A45 (7 ) # 
A55(7), 

A65 (7 ) * 

A75 (7 ) » 

A 6 5 1 7 ) * 

A 45 ( 7 ) # 

A105(7)» 
A 1 1 5 ( 7 ) 


C 

c 

c 


VlATFR (3t‘F* :T <400F ) 


t3 O 

Ms 

3 2 
o r? 

o £ 

r: r? 

> £ 

r- -r 

5rr 


000005 

DATA All( 1 ) * 

421(11/ 



000005 

“* *' ~ 491.7# 

16*016/ 



000005 

0ATA A3K1J* 

A32 ( 1 ) # A33 ( 1 )# 

A3Mi># 

4 35(3 ) / 

000005 

* 1 • 3355636# 

-4,9575766-5# 0.# 

0 • # 

0*/ 

C000C5 

PAT A A 4 1( 1 ) # 

A4 2 ( 1 ) , A43(l)» 

MM1)/ 

AA5<1>/ 

000005 

♦ I?0O.550b# 

-5. 70551 5^-2,-4.454586 6-4, 

0 • p 

0 • / 

000005 

DATA A c 1 ( 1 ) * 

452(1 ), A 53 ( 1 ) , 

Ai> V(l)# 

A$5(l ) / 

000005 

* 14.199322# 

-6.5267262, -.81013069, 

0# t 

0./ 

000005 

DATA Aft 1 ( 1 ) » 

A6?(l), A63 ( 1 ),. 

n< 

465(1)/ 

000005 

* 54. 4*1766# 

2, 5 b6 2966-2,-3 .54721 26- 5, 

0*9 

C./ 

CC0005 ~ 

" DATA ’ 471(1), 

A 7 2 ( 1 ) , A73 ( 1 ), 

A7AU)» 

A 7 U ( 1 ) / 

000005 

* 7,4432132# 

-6.C175c47, -.799409*2, 

J • f 

o./ 

000005 

DATA. A P 1 ( 1 1 * 

A 6 2 ( 1 ) « A a 3 < 1 ), 

A*U 1 >» 

Af 5(1)/ 

000005 

* 52.475785# 

-.26276099, 5 .0332 7C t -4, -4 

.4U*23£-7# 

o/ 

000005 

DATA 401(11, 

A9£(l), A93 ( l ) , 

A9M n> 

A.«<1J / 

000005 

* -ID. 46466* 

1 .104256?, 0 « , 

0 • 9 

C./ 

C0CO05 

DATA AIDitl)# 

A 1 02 ( 1 ) » A1C.3U ), 

AUMDi 

A 1 C j ( 3 ) / 

C00005 

♦ -1.D545655* 

5.3j 2?94?‘ -3,-6. 44tC -<.£-6,2, 


<*./ 

000005 

DATA AlUll)* 

411211), Al 13 ( 1 ), 

aiim i )* 

A 1 1 5 ( 1 )/ 


0 
cn 

1 


26263-6026-RU-00 



/ 


RliNX COMPILER 

<vrp .2.3*1 

09/26/76. 13. C5 

. 23 . 

PRCf S 

CC0OO5 

C 

* -9.A37’AoE-3, 

9 . 71 7223( -5 , 

-2.23C757E-7, 2.117J45C-1C#- 

7.53CUE-14/ 


‘ C 

f 

AHMHNIA ( 

-107.9F<T<190F» 



000005 

i> 

DATA AU(2>» 

A2K2)/ 




00C005 


* ’5i.e» 

17.032/ 




000005 


DATA ATjL(2), 

432(2)# 

433(2)# 

434(2)# 

425(2)/ 

QGG005 


♦ 1.31# 

0.# 

0.# 

0.# 

C./ 

C00005 


DATA A A 1 ( 2 ) , 

A A 2 ( 2 ) , 

A43 i 2 i> 

444(2)# 

A45(?) / 

000005 


* 1.093?*iE+3, 

-2.4e2955E+0, 

4.976430E-3, 

-4 .4749676-6# 

0./ 

000005 


DATA A5H2)# 

A52 ( 2 ) # 

A53 < 2 ) # 

A S A ( 2 ) , 

A 55 ( 2 ) / 

000005 


* i .i9?7?4E4i# 

-A,92179et+0» 

2.065C16E-1, 

-7 . s 79E 97r r 2 * 

(• . / 

OOCC05 


DATA A*l(2)# 

Ae2(?)# 

463(2)# 

A’6.4 ( 2 ) , 

Ab5<?)/ 

CCC005 


* 7.DA37*6fc4l, 

**1 • 1 72405E -1 # 

1 .93 1 707 1-4# 

-1.6449135-7# 

^ C * / 

000005 


DATA A 71(2) f 

472(2)# 

A73(2 ), 

A74< 2) , 

A75 1 2) / 

00000 5 


♦ l.?660*?6£4l, 

-1.113379(41# 

2.993126E+0# 

-4.6*97696-1# 

o. / 

000005 


DATA A81U), 

A 8 2 < 2 ) , 

A 8 3 f 2 > , 

A 8 4 1 2 ) » 

4H5(2) / 

000005 


* 3.5370AbEU» 

-2.996929E-1, 

6.623156t-4, 

-7,441f U5E-7, 3 

.5521546-10/ 

0QCQG5 


DATA A9 1 ( 2 ) , 

492(2)# 

A 93 ( 2 ) # 

494(2)# 

445(2)/ 

GO 00 05 


* .O’O ,e bE ♦ 3 » 

1.96609'>fc + 3» 

-4 .72671 5 E ♦ 2# 

5.054Cb6t ♦!#- 

2.v243t9e*0/ 

OOOCD5 


DATA At 01(2) , 

4102(2)# 

4103(2)# 

A10H(2), 

41 C5 ( £ ) / 

000005 


* -9.1601*66-1, 

3.999710E-3, 

-6.537242E-6, 

3 .04)94 35F-9, 

C.V 

000005 


DATA At 1 1 ( 2 ) , 

A 112(2)# 

All 3 ( 2 )# 

4114(2)# 

4115(2)/ 

000005 


+ 6.A?6501t-3» 

— 7.oC4641l — b # 

-7 .6997 59b -9# 

6 *02 35 33 E - i? # 

0./ 


V 

c 

-~C~' 

WCTHYL ALCCMOL (-140MT<3£0T) 



000005 

DATA Al 1 ( 3 ) f 

423(3)/ 




000005 


♦ 322.7, 

32.042 / 




C00005 


DATA A’l(3)» 

A3r (3># 

433(3)# 

434(3)# 

435 ( 3 )/ 

CC0005 


* t .203, 

C • » 

0 . t 

u. # 

( ./ 

GCC0C5 


PAT A AM(3U 

A 4 2 ( 3 ) # 

443(3)# 

444( 3) * 

4 4 6 (?) / 

000005 


* A . 7 OOP A bE + 2 » 

-2.4731C65+G# 

6.416fc29k-3# 

-7 • 0041 95h -6# 

?. 214439fc* 9/ 

GCC005 


DATA A5U3), 

452 (3) # 

Aa3<3>, 

454(3)# 

45S (3) / 

000005 


4 1.5054U£*1# 

-9.24C6 30b+v» 

3 .3661 36E+G# 

-1. 969200- *C # 

3. 3 69t 5 r f -.1 / 

000005 


DATA A *> 1 ( 3 )' , 

A 1 2 C 3 ) * 

A63 m» 

464(3)# 

46:. 13 > / 

0GCQO5 


* t.^r’ntjtfi, 

2.«J2A36 .-1 . 

-f .417f*72t-4# 

9 .?( 1 2 e n-7#- 

4. jCU'< c*1 3/ 

OOCOC5 


nATA A?l(3), 

472(3)# 

473(3)# 

474(3)# 

/ 7 5 ( 3 ) / 


o 

o> 
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09/28/76. 13. 05. 23. 


005 

cot 

•005 

005 

•005 

* 1 . 593164 E+ 1 , 

DATA 481 ( 3 ), 

"‘*"•..2 7 2853636 + 1 , • 
DATA 491 ( 3 ), 

* 1 « 9225966 + 2 , 

>00 5 

data 

4101 ( 3 ) , 

i 005 

♦ Q.Q 44433 l-Z$ 

>005 

DATA 

41 ) 1 ( 3 ), 

1005 

* 5 . 790525 E- 3 ,' 

C 

... 


C 


21 

c 



•005 

oata 

411 ( 4 ), 

)005 

♦ 

? 48 . 7 , 

1005 

DATA 

431 ( 4 ), 

>005 

♦ 

1 . 175 , 

>005 

OATA 

A 41 (Mi 

)005 

* o. 687825 E*l, 

)005 

OATA 

451 ( 4 ), 

>005 

♦ 3 . 270732 E+C, 

3005 

OATA 

. 461 ( 4 ), 

3005 ‘ 

* 1 . 332756 E+ 2 , 

3005 

DATA 

47114 ), 

0005 

♦ ®»$ 3 A 3 ? 2 E 4 l, 

3005 

TATA 

AP 1 ( 4 ) , 

0005 

+ -P. 34747 «, 6 * 0 , 

0005 ' 

DATA 

A 9 K 4 ), 

0005 

+ - 1 .P 385 P 8 E+ 3 , 

0005 

OATA 

4101 ( 4 ), 

0005 


0005 

data 

4111 ( 4 ), 

0005 

* f 

* - 5 . 24 P 9716 - 3 , 

c 

c 


ETHANE t- 

c 

000 £ 

’ " HAT A 

A U ( 5 ) ', 

0005 

* 

1 6 1 « £ s 

CG 05 

DATA 

4 * 1 ( 5 ), 


-2.109C9EE+1, 
AP2 (3 ), 
-1.153169E-1, 
492 ( 31 # 
-1.266769E+2, 
4102(3), 
2.0974176-4, 
4212(3), 
-i.404494E-5, 

C-55F<T<305F> 

A2U4) / 
102.93/ 
432(9), 
0 ., 
492(9), 

4 . 6 3 6 5 5 8 h - 1 , • 
452(4), 
1.573170b*),- 
402(4), 
-3.261757E-1, 
472(4), 
-T.6625751+2, 
462 (4 ), 
8. 5301 16E-2,* 
492(4)., 

1 . 1993666*3,- 
4102(4), 
-2.4035481-3, 
4112(4), 
4 .964 6 69 F-5 » ■ 


421(5)/ 
3C • 07/ 
432(5) . 


1.1443266+1,-4 

483(3), 

2. 3037956-4,-2 
493(3), 

3.1133446*1,-3 
4103(3), 
•6.31Q697E-7, 7. 

A113 ( 3 ), 
1.2056206-8, 3. 


.278643E+0, 
464(3), 
.155127E-7, 
494(3), 
.318410c +0, 
4104(3), 
394364E-10,- 
4)14(3), 
5166296-12, • 


433(4) 

0. 

A43 ( 4 ) 
1.6316856-3 
453(4) 
1.6079596*1 
' 463(4) 

1.111655E-3 
473(4) 
1 . 2526612*2 
483(4) 
2 .757696E-4 
493(4) 
2.944711E+2 
4103(4) 
5.7135126-6 
4113(4) 
1.1337476-7 


434(4) 

0. 

444(4) 
2.056547E-6 
454(4) 
5.2592436+0 
464(4) 
-1.611726E-6 
474(4) 
—4 .2656 6 2i+l 
484(4) 
3.6437246-7 
494(4 ) 
3.2150766+1 
4104(4) 
-6.3913021-9 
4114(4) 
9.2 356586-11 


433(5). 


434(5). 


PROFS . ' 

5.EC9906E-1/ 

465(3)/ 

7.551726-11/ . 

495(3)/ 

1. 325C51E-1/ 

4105(3 )/ 

-3.196966-13/ 

4115(3)/ : 

-8.670296-15/ 


A35(4>/ 

0 ./ 

A45 ( 4 ) / 

•1.0164466-9/ 

455(4)/ 

-6.2095016-1/ 

465(4)/ 

>.906674E-13/ 

475(4)/ __ 

5.3754636+0/ 

465(4)/ 

■1. 753176-10/ 

495(4)/ 

-1.3157266*0/ 

4105(4)/ 

2.65046c -12/ 

4115(4)/ 

-2 • 254 54 1 -1 4 / 


435 (5)/ 


o 


rffti-9209-69292 



>: • ., 


NX COH P H£R (VER.2. 3H) ... 09/26/76. 13.05.23,. . PROPS 


COO 5 * 1.16; O.f O.t 0., 0./ 

0005 DATA 4*115 ), A42(5), A43(5), A44(5I» A45(6)/ 

0005“ *~-T. 2760 f 4 E + 3, 4.5732546 + 1,- 1. 719 4 8 1 E- 1 , 2. 84G 43 9E-4, -1.7568696-7/ 

0005 DATA 451(5); AS2(5), A53(5), A540I, A55(5)/ 

0005 * 4, 5135206+1,-5. 8G3273E+1, 3. 388505E+1, -9. 165778E+0, 9.1547C4E-1/ . 

0005 ~~~6 aTA 461(5)/ A 1 2 ( 5 ) , A63(5), 464(5), At5(5)/ 

0005 * -3.433014E+2, 3 .901C41C+O-1 .478627E-2, 2 .45 1 166E -5, -1 , 5 1 86 1 2 6-6 / 

C005 HATA A7K5), A72(5)> A73(5>, A74(5), A75<5)/ 

[0005 ' * 9. 83108 OE+1,-1 .4637316+2, 8.4229286+1, -2. 191641E+1, 2.1298t3E+Q/ 

>0005 DATA 481(5), Ae2(5>, A83(5>» A64(5), A6515)/ 

J0005 *-1.7239436+1, 1 .9319206-1, -7. 9534226-4, i,385103E-6»-6.9C506E-10/ . 

10005* "DlTA" 401(5), 492(5), A93(5), 494(5), A45I5)/ 

10005 * 2.9998 *56 +4, -2 .01 74356+4, 5.C85813E+3,-5.697C99t+2» 2.342c?0£ + l/ 

J0005 DATA A J01 ( 5 ) , A1C215), 4103(5')# A1C4J5), 4105(5)/ 

10005’ ' *~-l .1428606 +0, 1 . 3170966-2,-5. C72525E-5, 8.39029<[e*-8*-5.1ieeGt-] 1/ 

10005 DATA Alll<5) » A112I5), A113(5)» 4114(5), 4115(5)/ 

30005 ♦ 1.123709E-2,-e. 3396226-5, 2. 759121E-7,-4. 393436-10, 2.65146+6-13/ 

"c 

c ixPTMANt (-250F<T<-12CF) 

00005 DATA* 411(6), ' AZ1 (6) / . 

00005 * 163.2, 16.04/ 

00005 DATA 431(6), 432(6), A3, %>, A34(6), A35(6)/ ... 

• 00005* * 1.32, 0., 9., 0., 0./ 

00005 OAT A A 4 1 ( 6 ) , 442(6), A43(6), 444(6), A45(6>/ 

00005 * -1.1?40O1E*3, 2. 4251425 +1,-1 .5893076-1, 4.547110e-4, -4. 9067146-7/ 

00005"" DATA "451(6), ’ A52(6), A53(6), A54(6>, A55(6)/ 

00005 * 1.365f»46+0, 8. 5576176+0, -3.74657CE+0, 5.8032476-1,-3. 2571586-2/ 

00005 CATA 461(6), 462(6), A63(6>, A64(6), A65(6)/ 

00005 * 'i;4AO««5iE+l, 3/bl4831‘-l, -3. 2235426-3, 1 ,'J 76420* -* , -1 . 35 31 2 36-3 / 

C00C5 DATA 471(6), A72(6), A73(6), 474(6), 475(6)/ 

C0C05 _ + 6.38ir‘»2t + l,-5.445C63t + l, 1 .8374936 + 1, -2.749C46 1 + 0, 1.5675456-1/ 

00005 'DATA 481(6), 482(6), ' A83I6), 4fc4(b), Ac6(6)/ 

00005 * -9. 5264P6E +0, 2 »U24 132€-1,-1 .5406206-3, 4.719715E-6, -5. 2116756-9/ 

00005 DATA A9 4 ( 6 ) , A92(6), A93(6), 494(6), A45(/)/ 

00005 " ' * 6, 1 1 06 3 1 E + 3', -6 • 76.6 5 2 9 E ♦ 2 , 1. 6805226+3,-2. 32iifcSc+2, 1.C744C9E + 1/ 

'CCbOS DATA 4)01(6), A102 ( 6 ) , Al03(o), A1C4(6), Alu5(6)/ 

*00005 * 3.48/ 4?eE-i,-l. 7209695-3, 3.t99297t-7, l.fc4C747 1 :-e,-4. 73323t-ll/ 
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C9/28/76. 13,05,23.. 


PROFS. 


00005 

C0005_ 


C 

C 

_s_ 


DATA All 1(6)# A112I6)# A113I6)# AllA(fc)# A115C6)/ 

*_ . , -5 .A639A7E-5, 2,ft562ACt-7,-7..fii.?C0E-10 Jf 8.lAtAS5E-13/ 

NITROGEN <-3A0F<T<-250F> 


00005 

•ocoos 

100005.. 

100005 

100005 

>00005.. 

>00005 

>00005 

>00005 

300005 

300005 

000005 

CGC005 

000005 

C00005_ 

000005 

000005 

000005 

000005 

000005 

00 COO 5 

000006 

000007 

C00011 

00001 3__ 

000015 

000017 

000021 

000023 ' 

000027 

CC0031 


DATA A1K7), 

* 113. 9# 
OATA A 31(7) » 

* l.AO# 

OATA AAII7)# 

* ?*<,A897AE*1#« 
’ "p A T A A 5 1 1 7 ) # 

* 3.R17173E+1#- 
"ATA A61(7>» 

“♦ ”7./0P716E*Ia‘ 
OATA A?1C7># 

* 2.102P02E+1#- 
~ OATA A81(7># 

* -1.139709E+1 , 
OATA 401(7)# 

~V 1. 71 66706 + <>» 
DATA Al01(7)# 

*■ 1.178000E-1#‘ 
OATA A 11 1 (7># 

* 1 .A3A0^1E-3» 


A21 ( 7} / 
28.016/ 
A32 ( 7 ) # 
0 .# 
AA2 ( 7 ) # 
■2.3G5556E-1*. 

A52 (7 ) . 
•1 « A 3128 91' ♦ l# 
A 6 2 ( 7 ) » 
■3.323232E-1# 
A72 ( 7 ) # 
•7.503727E40# 
Ab 2 1 7 ) # 
3.68G937E-1#- 
A92C7), 
■1.3719911 *A t 
A102I7)# 
-7.992A2AF-5#' 
A 112(7)# 
-3.768939F-6. 


A33 1 7 ) 

0. 

AA3(7> 
5...11.7599E-3 
A53 (7 ) 
3.C6A76AE+0 
A63 ( 7 ) 
2.281A69E-3 
a73(7J 
9.613273E-1 

A33 ( 7 ) 
■3.88G9A3E-3 
A93(7J 
A.1039A5E+3 
A 1 0 3 ( 7 ) 
■l_.A01615E-6 
A113 (7 ) 
•A. 379371 E-8 


A3 A ( 7 ) # 
0.» 
AAAI 7) » 
-2.3AG715E-5#, 
A 5 A I 7 > ,. 
-3.133777 t -l. 

A6A ( 7 ) # 
-7.A78632E-6# 
A7A ( 7 ) # 
-A.861U&E-2# 
A 6 A ( 7 ) » 
1.68701AE-5.- 
A9A< ?)» 
-5.A53515E+2# 
A10A ( 7 ) # 
0.# 
A11AI7I# 
1.2703966-10# 


I P(L. EO.O ) GC TO 20 

T2 • T*T 
T 3 " *■ T2+T 
TA » T2*T2 
Tb * 1900. /T 
TR? » T°*TR 
TP 3 » TR 2 +TR 
TOA ■ TP?*TR2 

At T» At. nr ( i ) 

AtT’«AlT*ALT 

ALTi.AtTR.ALT 


A35 ( 7 ) /. 
0 ./ 
AA5 17 > / 

. ..Oil. 
A55 ( 7 ) / 
1. 17t AA9E-2/ 
Ae5(7) / 
0 ./ 
A75I7)/ 
0 ./ 
Ae5<7>/ 
•2.6C92fcfcE“8/ 
A95 17) / 
2. 71662A E+1 / 
A105I7)/ 
0 ./ 
A1 15 1 7 ) / 
C./ 



* 

o , 
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V 

s 




*' 




jNX COHP Il E p tVg p »l*3M> . 


09/26/76. 13.05.2_3_*. 


PROPS. 


30033 


C 

C 


0 0035 


00037 

OOOA1 

00053_ 

00065 

00103 

00115 

C0133 

C0151 

C0167 

0626 r 

00213 

00215 

•00217 

10022 0_ 

JC0220 

>00221 


20 


AlT«»ALT?*AlT2 
" "fTu'IO PROPERTIES 

TF_» AI1 1U . - 

Ymw»A21(L> 

5H6V • A31(L) + A32(L)*T+A33Cl 1 AT2+A35 Ct)*T3*A35IL)*T6 > 

HFO » A41 (t)4AA2(L )*T + AA31L)*T2+AA^R)*T3+A<»5CL >+T.A . 

PV «T pV®~( A 5 1 1 U ♦ A 5 i U ) *TR4A53(t)*TR2+A5A(L)*TR3+A55(L)*TRA) 
PHOL • AAHl>*A62U>*T4A63<l>*T2+A6Ml>*T3MfcV<0+T5 

RHOV * PXP(A71(U+A72(L)*TR4A73UH1R24A7A(U*TR3+A75(U*TR5 

y I •" rvp(A61(L ) + AP2 C L 1 *T + A83 (L )*T2+A6A(l )*T3+A95(L)*TA) 

VI$V * C XP ( A9 j. (L ) +A9Z ( L ) *ALT + A93 ( L )* AL T2* A9A < 1 ) *AL T 3 + A95 ( L ) * 

YKL • A101(L>+A102tL>*T4A103(L>4T2+A10A(U+T34A105<U*T5 

- s r . All 1 (L )♦ A 112(1 ) *T + A113(L)*T2+AllA(L)4T34All5(C>*TA 
VISV * VISVM. 169750566 
vm»Vm/A.1697&0AEH 
RETURN 

rnwTTNUP ... ... . . -- •- 

RETURN '■ 

fNO 


> ..... 

ALTA) 


I 



o 
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ttjNX CO MPILER <Vgg«2.3H> . 09/28/76... 13.05. t 2.3». 


100016 
>00032 
>000 32_. 
>00032 


Situ o ruT T^e SANMA»P»wO»wl»Vi2#W3»VC^Vl»V2>y3»AA) 
p« «CPT f { V1~V3 ) **2+ ( Wl“W3 )**Z) 

TME^A* ACOS ( ( <W3“*«l>* CW2-W1 ) 4 ( V3-V1 > * < V2-V l > > / 

1 ?QRTtJ.U3-Vl)**2:MV3-Vl)**2). .. 

2 V((W2-K1 )**2+( V2-V1 ) **Z ) H 


>00067 
>00071 
>00107 
>00112 
>00130 
>00 i 3 A 


l*«R*TMFT A 

A0’ , »». 5* AbS < <W2*V3-V2*.w3)-(WC*V3-V0*h3 >♦ (WC*V2-VC*W? ) ) ... 


A12P3*.5*THETA*R*« 

A123-.5*ABS( ( Vi2*V3-V2*V»3)-(WlAV3-Vl*W3 


»(Ml*V2~Vl*w2> > 


A» A023+A1 2.B3-A12 3-AA ._ 
PETliPM 


>00135 W 






26263- 6026-RU- 00- 



ftU NX COMPI LER (VPR.gtlW) „ 09/29/76* 13.05_.23, 


SlinpnUTTM£ S1HPSN1 X, Y, AA, SF,AIPHA,DAIPHA, XC» YC»ALPKA0> 
000014' " CALL ~D(* V I ALPHA, 0X1, OS1»OA1,X#Y,XO,YO,ALPHAO) 

000032 ALRHA»ALRHA+DAIPHA/Z. 

000035 Xl«X+J0Yl/2.)+0ALPHA _ ... - — 

000041“ CALL DRV( ALPHA, QX2,DS2»DA2» XI, Y,XC»YO,ALPHAO> 

000057 ALPHA-ALPHA + CAlPHA/2 . 

C00062_ . x 2^X + 0X?*DALf H£ . ....... ... 

000060 CALL DRV( ALPHA, DX3, DS3,0A3, X2,Y» XO,YO»AL PHAO ) 

000103 X«X+f 0X1+4 • +DX2+DX3 ) + DAL PHA/6 • 

00011 2-. AA*AA+ ( 0 A 1+4* 40A2+0A3 >*0ALPHA/6* 

000121 ~ SF*SF+fn$i+4.+DS£+0$3)*0ALPHA/fc. 

C00130 R r TUPM 

0C 0131 P M 0 ... .... - - r*“ — 
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